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Chapter  1 


INTRODUCTION 

This  thesis  addresses  the  stabiUty  of  coastal  rubble  mound  stone  armor 
layers  exposed  to  water  waves.  The  most  common  coastal  mbble  mounds  are 
breakwaters,  jetties,  and  revetments.  Breakwaters  usually  have  a  large  part  of  their 
length  oriented  perpendicular  to  the  direction  of  wave  travel  and  their  primary  purpose 
is  to  produce  a  quiescent  area  in  their  lee  for  ship  anchorage  or  beach  sheltering.  Jetties 
are  usually  oriented  parallel  to  the  direction  of  wave  travel  and  are  used  to  produce  an 
artificial  channel  or  inlet.  This  channel  is  usually  an  entrance  to  a  harbor  and  is  used  for 
commercial  or  leisure  water-craft  navigation.  Revetments  typically  armor  coastal  or 
riverside  slopes  that  would  otherwise  erode  when  exposed  to  water  waves.  These  waves 
could  be  a  result  of  local  or  distant  storms  or  they  could  be  due  to  ship  wakes. 

Breakwater  failure  can  occur  due  to  a  number  of  different  failure  modes. 

The  dominant  failure  modes  are  shown  in  Figure  1.  Of  these,  seaside  armor  stability  on 
a  traditional  rubble  mound  is  critical  to  the  integrity  and  functionality  of  breakwaters 
and  is  tha«fore  the  focus  of  the  present  study.  Breakwater  stone-armor  layer  stability  is 
unique  from  other  structural  design  in  that  it  is  highly  variable,  and  this  variability  is 
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Figure  1.1.  Breakwater  failure  modes 


complex  and  difficult  to  quantify  for  the  wide  range  and  many  permutations  of  the  many 
variables  involved.  The  variability  is  dependent  on  the  stochastic  nature  of  both  the 
araior  boundary  conditions  and  loading.  Boundary  conditions  are  the  points  of  contact 
underneath  and  between  the  armor  stones.  The  boundary  conditions  are  uncertain 
because  of  the  irregular  stone  shape  and  stone  placement.  The  geotechnical 
characteristics  of  the  foundation  are  also  often  uncertain.  The  loading  of  engineering 
interest  is  primarily  due  to  incident  storm  waves  that  vary  with  storm  intensity  and 
storm  location.  The  impact  of  the  waves  varies  with  water  level  and  local  bathymetry. 
The  local  water  level  varies  with  tide  and  storm  surge.  Waves  can  dislodge  armor  units 
by  uplifting,  rolling,  or  sliding  individual  units  or  by  causing  en  masse  movement  of  the 
entire  armor  layer.  Armor  units  can  be  dislodged  from  the  upslope  layer  without 
jeopardizing  the  integrity  of  the  armor  layer.  But  if  enough  armor  units  move  and  the 
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underlayer  is  exposed,  then  the  breakwater  can  erode  quickly.  This  structural  integrity 
threshold  is  uncertain  because  damage,  defined  by  the  eroded  volume,  may  be  focused 
at  a  point  or  distributed  over  a  broad  area.  Because  the  toe  forms  a  foundation  for  the 
armor  layer,  if  the  toe  armor  is  mobile,  then  the  entire  armor  layer  can  mobile.  How  the 
armor  layer  responds  to  armor  movement  at  the  toe  is  quite  variable. 

The  variability  in  both  the  loading  and  boundary  conditions  demands  the 
consideration  of  the  randomness  of  each  stochastic  parameter.  Many  clients,  including 
the  U.S.  Army  Corps  of  Engineers,  now  mandate  that  a  risk  analysis  be  done  as  part  of 
major  breakwater  rehabilitation  studies.  This  requirement  is  intended  to  provide  a 
standardized  engineering  economic  analysis  technique  for  comparing  competing 
alternatives.  The  engineering  performance  study  that  must  form  the  basis  of  the  risk 
analysis  is  often  accomplished  using  a  reliability  analysis,  where  the  reliability,  or 
conversely,  the  probability  of  failure,  is  quantified  for  each  alternative.  In  the  reliability 
analysis,  all  alternatives  that  fall  below  a  predefined  level  of  reliability  are  either 
modified  or  discarded.  For  example,  a  breakwater  may  have  several  alternative  armor 
layer  designs  including  randomly-placed  stone,  pattern-placed  stone,  and  randomly- 
placed  concrete  armor  units.  Each  of  these  alternatives  will  have  a  different  associated 
cost  and  a  different  probability  of  failure  for  the  design  single  storm  or  sequence  of 
storms.  Determining  the  probability  of  failure  through  instability  of  the  armor  layer 
requires  knowledge  of  the  rate  of  deterioration  or  the  rate  at  which  stones  are  displaced, 
which  has  not  been  quantified  for  breakwater  armor  layers.  Reliability  methods  have 
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only  recently  been  adapted  to  breakwater  design,  and  very  little  data  exist  to  support 
prediction  of  rubble-mound  deterioration.  Despite  the  many  hundreds,  and  perhaps 
thousands,  of  studies  on  breakwater  armor  stability,  there  have  been  few  generalized 
studies  of  long-term  deterioration  of  traditional  mbble-mound  armor  layers.  Moreover, 
there  have  been  no  generalized  studies  of  deterioration  due  to  variations  in  storm 
sequences  using  random  waves.  Finally,  damage  experiments  to  date  have  been 
primarily  conducted  with  nonbreaking  waves,  which  is  atypical  of  many  environments 
in  which  breakwaters  are  constructed. 

The  purpose  of  this  study  was  to  first  identify  the  primary  mechanisms  of 
stone  instability  and  damage  development  when  the  breakwater  armor  layer  is  exposed 
to  depth  limited,  breaking  waves.  Chapter  2  discusses  initiation  of  stone  movement.  A 
relation  for  predicting  the  initiation  of  stone  movement  when  exposed  to  vertical  uplift 
of  normally-incident  plunging  breaking  waves  is  presented  in  Chapter  2.  The  equation 
is  verified  using  data  from  a  small-scale  two-dimensional  flume  physical  model  study. 
This  portion  of  the  research  effort  qualitatively  addresses  incipient  motion  of  stone 
armor.  This  section  includes  the  traditional  development  of  a  stability  prediction 
equation  and  some  insight  into  the  effect  of  wave  steepness  on  instability. 

A  primary  goal  of  this  study  was  to  establish  predictive  relations  for  damage 
development  on  traditional  breakwater  sections  for  single  storms  and  for  storm 
sequences  given  depth-limited  normally-incident  waves.  Chapter  3  introduces  the 
subject  with  a  discussion  of  historical  stability  and  damage  development  physical  model 
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measurements  and  the  relations  derived  from  these  past  studies.  In  Chapter  4,  a  new 
flume  study  is  discussed  where  breakwater  profiles  were  measured  as  damage 
progressed  on  a  breakwater  cross  section  exposed  to  normally-incident  depth-limited 
breaking  waves.  A  new  device  for  measuring  breakwater  profiles  is  discussed  in  some 
detail.  New  parameters  are  defined  for  prescribing  the  engineering  characteristics  of  the 
eroded  profile.  In  Chapter  5,  measurements  from  this  study  are  discussed  in  detail.  In 
Chapter  6,  profile  and  wave  data  are  analyzed  to  produce  spatial  and  temporal  relations 
for  predicting  the  mean  and  standard  deviation  of  eroded  area,  eroded  depth,  eroded 
length,  and  remaining  cover  depth  on  a  breakwater  that  is  exposed  to  normally-incident 
depth-limited  wave  conditions.  These  equations  should  be  valuable  in  support  of 
reliability  analyses  as  part  of  comprehensive  risk  analyses. 
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Chapter  2 


INCIPIENT  STABILITY  OF  BREAKWATER  ARMOR  UNITS 


2.1  Armor  Stability  Equations 

Extensive  research  on  breakwater  armor  stability  has  produced  many 
empirical  stability  models.  PIANC(1976)  provides  a  summary  of  early  stability 
models  based  on  regular  wave  experiments.  PIANC  (1992)  provides  a  discussion  of 
more  recent  irregular-wave-based  models.  The  most  widely  known  empirical  stability 
model  was  developed  by  Hudson  (1958, 1959),  following  the  pioneering  work  of 
Iribarren  (1938)  and  is  typically  seen  in  the  following  form 


W  = 


K^iS^- D^cotd 


(2.1) 


where 


W  =  weight  of  armor  unit 

Yfl  =  specific  weight  of  armor  unit  material 

H  =  design  wave  height  at  stracture  toe 

Kj)  =  tabulated  empirical  stability  coefficient 

5^  =  specific  gravity  of  armor  unit  material 

0  =  seaside  angle  of  armor  slope  relative  to  horizontal 
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Hudson  also  expressed  this  equation  in  a  slightly  different  form  as 


=  (/ToCOte)*^  = 


H 


(2.2) 


where  is  the  dimensionless  stability  number  and  is  termed  the 

nominal  stone  diameter.  Hudson  did  not  use  the  nominal  stone  diameter  variable.  It 
was  used  by  both  Thomson  and  Shuttler  (1976)  and  Van  der  Meer  (1988)  who  noted 
that  Z)„5o  is  the  length  of  a  side  of  a  cube  with  volume  equivalent  to  that  of  the  median  in 
the  stone  weight  distribution.  A  glance  at  equations  2.1  and  2.2  shows  that  stone 
movement  due  to  wave  forcing,  characterized  by  the  wave  height,  is  resisted  primarily 
by  the  stone  weight.  As  shown  by  Hudson  (1958),  Ahrens  and  McCartney  (1975),  Van 
der  Meer  (1988)  and  others,  the  wave  shape,  structure  porosity,  armor  shape,  and 
stmcture  cross-sectional  shape  all  affect  the  wave  force.  Hudson  also  noted  that  the 
resistance  to  movement  is  affected  by  the  fiiction  between  armor  units,  stone  shape, 
upslope  armor  layer  weight,  and  armor  slope.  Equation  2.1  is  preferable  from  a  physical 
perspective  because  it  maintains  the  inertial  form  of  the  forces,  which  is  appropriate  for 
armor  stability.  But  equation  2.2  combines  the  stability  coefficient  with  the  stracture 
slope,  which  is  desirable  because  the  stability  coefficient  is  a  function  of  the  structure 
slope  (e.g.  Shore  Protection  Manual  (SPM)  1984).  This  second  form  also  presents  the 
stone  stability  equation  in  a  form  that  is  similar  to  sediment  transport  formulae.  This 
will  be  discussed  further  later  in  this  chapter. 
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In  equation  2.2,  the  stability  number  is  a  convenient  scale  for  characterizing 
the  incident  wave  height  relative  to  armor  stone  size.  It  is  typically  in  the  range  of  1  to  4 
for  stable  or  marginally  stable  rubble  mound  breakwaters  (Van  der  Meer  1988).  The 
SPM  (1984)  provides  values  for  various  armor  layers  at  different  slopes  on 
breakwater  trunks  and  heads  exposed  to  breaking  and  nonbreaking  waves.  Zero-damage 
Kq  values  are  typically  used  corresponding  to  less  than  2  percent  by  count,  or  five 
percent  by  volume,  of  the  armor  in  a  layer  being  displaced.  Using  this  no-damage 
guidance  from  the  SPM  (1984),  the  stability  number  is  in  the  range  of  1  to  1.6  for 
angular  stone  armor  layers  at  slopes  of  1  V:2H  or  steeper  exposed  to  breaking  waves. 

The  SPM  specifies  values  up  to  2.2  for  stone  armor  layers  ait  slopes  of  1  V:3H 
exposed  to  nonbreaking  waves.  So  for  most  stable  coastal  breakwaters,  the  stability 
number  covers  a  narrow  range  from  1  to  2.5.  For  deformable  stractures  where  the 
armor  stone  is  expected  to  be  mobile,  such  as  S-shaped  breakwaters  and  berm 
breakwaters.  Van  der  Meer  (1988)  suggests  N^  =  3  -6. 

Hundreds  of  studies  have  been  conducted  to  quantify  the  single  empirical 
parameter  Kp  in  equation  2. 1  for  the  wide  vaiiety  of  prototype  conditions  that  might 
exist.  Originally  Hudson  only  explicitly  included  the  effect  of  regular  wave  height, 
structure  slope,  and  armor  stone  specific  weight.  Hudson  found  no  clear  effect  of  wave 
period  on.armor  stability  for  the  nonbreaking  regular  wave  conditions  he  and  his 
colleagues  tested.  They  simply  determined  Kp  corresponding  to  the  lowest  stability 
condition  over  a  range  of  typical  wave  periods.  The  Hudson  equation  has  been 
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extended  to  include  the  effects  of  irregular  breaking  and  nonbreaking  waves  and  wave 
period  (Ahrens  1975,  Ahrens  and  McCartney  1975,  Carver  and  Wright  1991,  SPM 
1984),  various  armor  layer  types  and  armor  gradation  (e.g.  SPM  1984),  and  number  of 
waves  (Medina  and  McDougal  1988). 

Other  regular-wave  stability  formulations  that  have  been  utilized  in  recent 
years  include  Hedar  (1960, 1986),  Ahrens  and  McCartney  (1975),  and  Losada  and 
Gimemenz-Curto  (1979).  Ahrens  (1975),  Ahrens  and  McCartney  (1975),  Losada  and 
Gimemenz-Curto  (1979),  and  Pilarczyk  and  Den  Boer  (1983)  all  showed  dependence  of 
wave  period  on  stability  number  for  regular  waves.  Each  of  these  authors  showed  that 
minimum  stability  occurred  for  surf  similarity  numbers  or  Mbarren  nmnbers  between  2 
and  4,  corresponding  to  plunging  to  collapsing  breakers.  The  surf  similarity  number 
was  defined  as  ^  =  tan9  /  {H/Lf  where  6  =  structure  slope,  H  =  regular  wave  height, 
and  L  -  local  or  deep  water  wave  length.  Thompson  and  Shuttler  (1975)  conducted  an 
extensive  series  of  irregular-wave  armor  stability  experiments.  Their  conclusions  on 
damage  progression  were  insightful  and  are  discussed  in  the  next  chapter.  Using 
Thompson  and  Shuttler’s  data  and  data  from  his  own  experiments.  Van  der  Meer  (1988) 
developed  a  stability  model  which  explicitly  included  wave  period,  stmcture 
permeability,  storm  duration  and  damage  for  a  single  design  storm.  The  Hudson  and 
Van  der  Meer  equations  will  be  discussed  further  in  the  following  chapters.  The  above 
mentioned  stability  models  predict  minor  damage  reasonably  well,  although  poor 
predictions  are  common.  Pfeiffer  (1991)  compared  the  models  of  Hudson  (1958), 
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Hedar  (1986),  Losada  and  Gimemenz-Curto  (1979),  and  Van  der  Meer  (1988).  Pfeiffer 
found  that  the  Hudson  equation  and  the  Losada  and  Gimemenz-Curto  equations 
performed  the  best  when  compared  to  field  data,  but  none  of  the  equations  matched  the 
prototype  records  all  that  well.  Pfeiffer  also  noted  that  none  of  the  stability  models 
mentioned  above  are  specifically  suited  for  predicting  extended  damage. 

2.2  Armor  Incipient  Motion  Studies 

The  empirical  stability  models  of  Iribarren  (1938),  Hudson  (1958),  and 
others  are  based  on  a  free  body  analysis  of  an  armor  unit  undergoing  forcing  due  to 
shallow- water  waves.  Early  stability  models  assumed  1)  the  principal  wave  force  was 
due  to  down-  or  up-rush  on  an  unsheltered  and  unrestrained  imit,  2)  the  wave  force  was 
drag  dominant,  and  3)  the  drag  force  would  be  critical  if  the  maximum  horizontal  fluid 
velocity  was  used,  which  was  considered  to  be  proportional  to  the  shallow-water 
incident  wave  celerity.  But  for  an  intact  structure  and  prior  to  initiation  of  incipient 
motion,  the  armor  units  are  typically  partially  hidden  and  restrained  from  up  or  down 
slope  movement;  so  lift,  inertia,  and  convection  across  the  armor  layer  must  be 
considered.  Moreover,  Sawaragi  et  al.  (1982)  showed  that  the  maximum  fluid  velocity 
on  a  rubble  mound  was  not  necessarily  proportional  to  the  wave  celerity.  Sigurdsson 
(1962)  made  force  measurements  on  armor  composed  of  spheres  set  on  extremely  steep 
slopes  with  an  impermeable  underlayer  and  derived  incipient  equations  of  motion;  but 
concluded  by  stating  that  the  dominant  mechanism  of  initiation  of  armor  motion  was 
still  unknown  and  required  further  investigation.  Although  many  authors  have 
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measured  wave  kinematics  on  armor  layers  and  the  resulting  forces,  including  Mizutani 
et  al.  (1992),  Torum  and  Van  Gent  (1992),  Torum  (1994),  and  Cornett  and  Mansard 
(1994),  there  have  been  few  observations  of  incipient  movement  of  armor  units 
discussed  in  the  literature.  As  such,  the  relationships  between  incipient  motion,  wave 
kinematics,  and  forces  on  armor  units  are  still  unknown. 

Melby  (1987)  and  McDougal  et  al.  (1988)  discussed  a  model  for  predicting 
the  wave  forces  on  dolos  concrete  armor  units  and  the  resulting  incipient  motion  of  a 
lone  dolos  in  one  of  two  orientations.  Their  model  utilized  linear  wave  theory  and 
Morison  forcing  (Morison  et  al.  1950)  with  added  mass  coefficients  for  wave  slamming 
from  Kaplan  and  Silbert  (1976)  and  Kaplan  (1979).  Kobayashi  et  al.  (1990)  presented  a 
numerical  model  for  predicting  the  displacement  of  armor  on  a  traditional  rubble 
mound.  The  shallow  water  wave  model  interacted  with  a  permeable  flow  model  and 
hydrodynamic  drag,  inertia,  and  lift  forces  were  computed  using  a  Morison  formulation. 
The  model  was  limited  to  forces  parallel  to  the  sfructure  because  only  depth-averaged 
velocities  were  predicted  by  this  one-dimensional  model.  Torum  and  Van  Gent  (1992) 
discussed  a  similar  wave  model  and  compared  it  to  velocity  measurements  above  a  berm 
breakwater.  Torum  (1994)  discussed  the  measurements  further.  Although  two- 
dimensional  velocities  were  measured,  vertical  flow  in  the  breaking  wave  was  not 
modeled  numerically.  In  addition,  Torum  noted  that  the  inertial  force  was  not  well 
defined  by  the  traditional  inertia  term  of  the  Morison  equation.  Cornett  and  Mansard 
(1994)  described  an  experiment  where  forces  were  measured  on  a  panel  of  stones.  This 
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approach  was  unique  and  yielded  insight  into  the  average  frictional  force  on  sections  of 
the  armor  layer.  They  found  that  the  frictional  force  tending  to  dislodge  armor  units  was 
greatest  below  the  still  water  level  and  that  the  character  of  the  force  depended  strongly 
on  the  type  of  wave  breaking. 

This  chapter  discusses  a  series  of  physical  model  experiments  to  identify 
and  develop  predictive  models  for  breakwater  armor  incipient  motion  and  to  relate  this 
motion  to  existing  empirical  stability  relationships.  The  experiments  were  conducted  in 
wave  flumes  at  the  Waterways  Experiment  Station  in  Vicksburg,  Mississippi.  The  first 
experiment  consisted  of  measuring  wave-induced  fluid  velocities  on  and  within  the 
armor  layer  and  runup/down.  In  addition,  free  surface  oscillations  were  measured  while 
observing  armor  motion  on  stone  and  Core-Loc  armor.  The  observations  from  this  early 
study  led  to  an  incipient  motion  experiment  using  a  fixed-sphere  armor  layer  with 
several  loose  spheres  placed  at  various  depths  within  the  armor  layer.  A  dominant 
incipient  armor  motion  mode  and  predictive  stability  equation  were  verified. 

2.3  Experimental  Setup 

The  initial  experiments  were  conducted  to  determine  the  nature  of  armor 
incipient  motion  and  surrounding  flow.  The  instrumentation  included  a  laser  Doppler 
velocimeter  (LDV),  high-resolution  video,  and  ranup  and  vertical  free-surface-piercing 
gages  near  and  within  the  armor  layer.  The  experiments  discussed  in  this  chapter  for 
incipient  motion  were  all  carried  out  using  regular  monochromatic  waves.  Data  analysis 
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was  performed  on  short  segments  of  between  five  and  ten  uniform  waves  to  develop 
clear  relationships  between  wave  parameters  and  armor  motion. 


The  primary  study  was  conducted  in  a  46-m  long  by  0.46-m  wide  by  1-m 
deep  flume,  with  an  offshore  slope  of  1V:30H  (Figure  2.1).  A  conventional  rubble 
mound  cross  section  was  constmcted  with  various  seaward  slopes  and  armor  types 
(Figure  2.2).  Table  2.1  lists  the  different  test  plans.  In  this  thesis,  only  the  sphere  and 
stone  armor  plans  will  be  discussed  (Test  Plans  3  and  4)  in  order  to  maintain  continuity 
with  the  stone  layer  damage  discussion  in  subsequent  chapters.  For  Plans  3  and  4, 
velocity  measurements,  sampled  at  100  hz,  were  made  throughout  the  water  column 
from  the  toe  to  two  armor  dimensions  above  the  still  water  level.  The  ranges  of  physical 
quantities  and  common  dimensionless  parameters  for  Plans  3  and  4  are  listed  in  Table 
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Sea  Side  Harbor  Side 


Figure  2.2.  Definition  sketch  for  typical  structure  profile 


Table  2.1 

Experimental  Plans  for  Incipient  Motion  Experiments 

Plan 

Armor 

Type 

Armor 

Weight 

W,g 

Nominal 

Armor 

Dia. 

D„=(V)’« 

cm 

Struc¬ 
ture 
Slope, 
cot  6 

Off- 
Shore 
Slope 
cot  a 

Runup 

down 

Meas. 

Force 

Meas. 

1 

Core-Loc 

220 

4.6 

1.33 

100 

no 

no 

2 

Core-Loc 

105 

3.6 

1.5 

20 

no 

no 

3 

Stone 

200 

4.6 

2.0 

30 

yes 

no 

4 

Sphere 

58 

3.03 

2.0 

30 

yes 

no 

5 

Sphere 

212 

5.6 

2.0 

30 

yes 

yes 
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Table  2.2 

Ranges  of  Measured  Physical  Quantities  and  Common  Dimensionless 
Parameters  for  Plans  3  and  4 

Parameter 

Range 

Wave  height  at  toe,  H, 

2.4  to  1 8  cm 

Wave  period,  T 

0.75  to  4.0  sec 

Water  depth  at  toe,  h, 

to  24  cm 

Wave  steepness,  5^=  H,/Lo 

Lg  =  deepwater  wave  length 

0.007  to  0.1 

Surf  similarity  parameter,  ^  =  tan  a/iH/Lg)^ 
tan  a  =  1/30  =  beach  slope 

0.1  to  1.4 

Surf  similarity  parameter,  ^  =  tan 
tan  0  =  1/2  =  structure  slope 

1.5  to  21 

Relative  depth,  h/Lg 

0.009  to  0.28 

Relative  wave  height,  H/h, 

0.06  to  0.88 

The  LDV  was  a  two-watt  argon-ion  two-component  device  assembled  by 
the  Dantec  Coiporation.  Here,  two-component  describes  the  fact  that  the  LDV 
measured  velocities  in  two  orthogonal  directions.  For  this  study,  these  directions  were 
always  vertically  upward  and  horizontal,  in  the  direction  of  wave  propagation.  The 
LDV  worked  in  the  back-scatter  mode  using  a  non-intrasive  probe  which  contained  both 
the  emission  and  receiving  optics.  The  benefits  of  this  device  included  nonintrasive 
measurements,  small  measurement  volume,  clean  drop-outs,  high  sampling  rate,  and  no 
required  calibration.  Drop-outs  were  situations  where  there  was  no  measurable 
backscattered  signal  (e.g.  when  the  laser  beams  crossed  at  a  point  above  water,  after  the 
wave  crest  had  passed  and  the  probe  was  above  the  level  of  the  wave  trough).  During  a 
drop-out,  the  LDV  would  produce  a  constant  signal  at  the  last  measured  value.  Because 
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only  the  peak  measured  velocity  values  were  used  herein,  drop-outs  did  not  pose  a 
particular  problem. 

For  this  experiment,  the  LDV  probe  was  pointed  through  the  glass  flume 
wall  into  the  voids  within  the  armor  layer  and  measurements  were  made  of  the  internal 
flow  within  the  porous  media.  Many  of  the  voids  between  armor  units  were  more  than  a 
nominal  armor  diameter  deep,  so  measurements  could  be  made  outside  the  flume  wall 
boundary  layer.  One  drawback  to  the  LDV  was  the  requirement  of  a  full  time  operator 
with  continuous  attention  to  detail.  Also,  because  of  the  small  measurement  volume, 
small  changes  in  measurement  location  in  the  sheltered  region  behind  a  stone  or  within 
the  armor  layer  often  yielded  large  variations  in  measured  average  peak  velocities. 
Therefore,  the  instrument  required  many  measurements  to  map  the  flow  field.  So  data 
analysis  requirements  were  substantial  for  this  experiment. 

The  wave  heights  were  determined  using  free  surface  measurements  from  a 
vertical  capacitance-type  gage  positioned  at  the  location  of  the  structure  toe  with  no 
structure  in  place.  Synthetic  rubberized  horse-hair  mats  were  placed  landward  of  the 
stmcture  location  to  absorb  the  waves.  The  sampling  rate  for  free  surface  measurements 
was  20  hz.  The  zero-downcrossing  wave  height  was  computed  as  the  average  height 
from  a  burst  of  approximately  ten  regular  waves. 
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2.4  Incipient  Motion  Observations 

Several  dominant  incipient  motion  modes  were  identified  during  the  stone 
stability  experiment  in  Plan  3,  Table  2.1.  The  following  descriptions  pertain  to  initial 
armor  motion  on  an  intact,  as-built  structure.  Rolling  was  the  only  mode  of  motion  for 
stones  on  the  toe.  Although  both  onshore  and  offshore  motion  was  observed,  the  toe 
units  always  moved  out  of  the  layer  in  the  onshore  direction.  Armor  near  the  still  water 
level  was  more  likely  to  displace  than  armor  in  other  areas.  This  appeared  to  be  due  to 
the  fact  that  the  armor  was  loosened  in  this  area  due  to  high  velocities  in  the  breaking 
wave  jet.  Once  loosened,  the  motion  would  depend  on  the  armor  shape  and  its  position. 
If  the  armor  shape  was  flat,  then  the  armor  unit  would  flop  back  and  forth  until  it  rolled 
out  of  the  armor  layer,  generally  rolling  upslope  during  uprush.  If  the  armor  shape  was 
rounded,  which  was  normally  the  case,  the  armor  units  would  jump  vertically  under  the 
steep  wave  face  if  the  wave  was  severely  plunging  or  collapsing.  If  the  wave  was 
surging,  then  loose  units  would  only  be  displaced  if  they  were  exposed.  There  did  not 
appear  to  be  sufficient  lift  in  downrush  or  uprash  flows  along  the  armor  layer  to  displace 
the  stones  unless  they  were  odd  shaped  (flat).  The  only  displacement  mechanism 
observed  for  rounded  stones  sufficiently  hidden  in  the  armor  layer  was  uplift  under  the 
steep  wave  face. 

These  observations  indicated  that  a  fluid  velocity  or  acceleration  component 
in  the  vertical  direction  is  normally  required  to  initiate  armor  motion  for  hidden  armor 
units.  Additionally,  this  early  qualitative  study  indicated  that,  for  a  given  wave  height. 
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incipient  motion  was  primarily  influenced  by  wave  shape,  stone  position,  stone 
exposure,  and  stone  looseness. 

2,5  Experimental  Measurements 

Throughout  the  experiments,  vertical  and  horizontal  wave  velocities  were 
measured  in  the  vicinity  of  the  armor  layer.  Figure  2.3  shows  typical  time  series  of  the 
horizontal  and  vertical  velocities  on  the  structure  measured  for  one  run  of  plan  3  with 
the  following  characteristics:  H,  =  12  cm,  T=  \s,d,  =  24  cm,  and  di  =  8.8  cm,  where 
is  the  average  wave  height  measured  at  the  toe,  fis  the  average  wave  period,  d,  is  the 
toe  water  depth,  and  ^4  is  the  depth  of  the  laser  measurement.  The  laser  measurement 
was  made  1  cm  outside  the  armor  layer,  measured  perpendicular  to  the  outer  armor  layer 
profile  line.  The  sign  convention  is  such  that  the  horizontal  velocity  is  positive  seaward 
while  the  vertical  velocity  is  positive  upward.  Typically  for  these  measurements,  the 
horizontal  velocity  signal  was  considerably  smoother  than  the  vertical  velocity  signal, 
due  primarily  to  the  small  amplitude  of  the  vertical  velocity  relative  to  the  horizontal 
velocity  and  the  relatively  large  amount  of  turbulence  near  the  armor  layer.  Figure  2.4 
shows  a  velocity  vector  time  series  over  one  wave  period,  measured  1  cm  outside  the 
armor  layer  profile  line  for  one  run  of  plan  3  with  the  following  characteristics:  H,  =  8.4 
cm,  r=  2  s,d,  =  24  cm,  and  d^=  17  cm.  Also  shown  is  the  wave  profile  at  the  point  of 
maximum  vertical  velocity.  The  plot  shows  a  large  vertical  velocity  just  below  the  steep 
wave  front.  Observed  maximum  stone  movement  for  this  wave  profile  position  is  also 
shown.  Figure  2.5  shows  an  example  of  vertical  velocities  outside  and  inside  the  armor 
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layer.  These  measurements  were  also  made  during  one  run  of  plan  3  with  the  following 
characteristics:  H,  =  12  cm,  f-  1  s,  4  =  24  cm,  and  4,  =  8.8  cm  for  the  external 
measurement  and  4  =  13  cm  for  the  internal  measurement.  The  measurement  locations 
for  these  time  series  are  shown  in  Figure  2.6,  where  the  stmcture  slope  is  1  V:2H.  Here 
it  is  clear  that  the  velocities  within  the  armor  layer  are  highly  irregular  due  to 
turbulence. 


Figure  2.3.  Velocity  time  series  for  one  wave  period  with  H,  -  12cm,  T=l  s,  d,  =  24  cm 
and  4  =  8.8  cm  measured  1  cm  outside  the  armor  layer 
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Figure  2.4.  VelcKity  vector  for  one  wave  period  for  armor  lifting  measured  1  cm  outside 
the  armor  layer  with  ti,  =  8.4  cm,  f-  2  s,  4  =  24  cm,  and  =  17  cm 
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Figure  2.5.  Vertical  velocity  time  series  outside  (top)  and  inside  (bottom)  armor  layer 
with  H,  =  12  cm,  f  =  1  s,  4  =  24  cm,  and  4  =  8.8  cm  for  the  outer  measurement  and  4. 
=  13  cm  for  the  iimer  measurement 
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Figure  2.6.  Measurement  locations  for  vertical  velocities  shown  in  Figure  2.5 

Figure  2.7  shows  a  typical  plot  of  vertical  velocity  measured  1 

cm  outside  the  armor  layer  at  various  depths,  versus  the  square  root  of  wave  steepness, 
where  Vj/j  is  the  average  of  the  highest  one-third  peak  velocities  for  the  burst  of  regular 
waves,  g  =  gravitational  acceleration,  structure  slope  =  1  V:2H,  and  =  deep  water 
wave  length.  For  simplicity,  v  is  used  instead  of  v,;3  in  Figure  2.7  and  hereafter. 
Relative  laser  depth,  rd  =  is  the  ratio  of  the  depth  of  the  laser  to  the  depth  at  the 
toe,  measured  from  the  still  water  level.  As  noted  by  Sawaragi  et  al.  (1982),  maximum 
non-dimensional  velocities  commonly  occurred  for  collapsing  to  plunging  breaking 
waves.  Pilarczyk  and  Den  Boer  (1983)  showed  minimum  stability  for  1V:2H  slope 
occurred  for  ^  =  3.3.  In  this  case,  this  point  occurs  at  =  0.15,  which  is  near  the 
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maximum  measured  velocity  for  the  waves  with  T=  2  s.  The  peak  vertical  velocities  for 


given  wave  period  decreased  with  increase  of  rd  in  this  figure. 


Vertical  Velocity  at  Armor  Surface 


°  T=1,  rd=0.36  -  T=2,  rd=0.36  -  T=1,  rd=0.5 
•  T=2,  rd=0.5  -  T=1,  rd=0.7  ▼  T=2,  rd=0.7 


Figure  2.7.  Maximum  vertical  velocity  versus  the  square  root  of  wave  steepness  for 
wave  periods  T=  1  and  2  s  and  relative  depths  rd  =  0.36, 0.5,  and  0.7 


2.6  Incipient  Motion  Prediction 

The  previous  experimental  results  indicated  that  one  of  the  dominant 
incipient  motion  modes  was  due  to  the  vertical  force  occurring  under  the  steep  wave 
front.  The  balance  of  forces  for  vertical  incipient  armor  motion  with  no  external 
restraints  yields  the  equality  between  the  submerged  armor  weight  and  the  vertical  fluid 
force;  W  =  F^.  The  vertical  force  at  the  steep  wave  front  can  be  described  by  the 
Morison  equation  (Morison  et  al.  1950). 
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(2.3) 


+  PVCJ^ 

where  p  =  fluid  density,  A  =  cross  sectional  area  of  armor  unit  in  direction  of  flow,  Q,  = 
drag  coefficient,  v  =  peak  vertical  velocity  just  above  the  armor  layer  as  plotted  in 
Figure  2.7,  V  =  armor  unit  volume,  C„  =  inertia  coefficient,  and  dv/dt  =  total  fluid 
acceleration.  The  drag  force,  the  first  term  on  the  right  side  of  Equation  2.3,  can  be 
expressed  as  a  function  of  the  armor  nominal  diameter,  D„,  by  introducing  an  armor 
shape  factor,  as  follows 

A  =  (2.4) 


(2.5) 


(2.6) 


At  the  point  of  maximum  vertical  fluid  velocity,  the  local  vertical  fluid 
acceleration,  dvidt,  and  horizontal  velocity,  u,  are  negligible.  As  such,  the  total 
acceleration  reduces  to  a  convective  term. 
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max 


(2.7) 


'  f  dv  dv 

-  =  -  +  U -  +  V - 


dv 

X  V - 

dy 


where  x  =  horizontal  coordinate  and  y  =  vertical  coordinate.  If  we  assume  the 
convective  acceleration  across  the  armor  layer  to  vary  linearly  vertically,  then  the 
acceleration  can  be  expressed  as 


dv 

dt 


max 


(2.8) 


where  is  an  empirical  coefficient  of  order  unity.  The  maximum  inertial  fluid  force  in 

Equation  2.3  can  thus  be  reduced  to 


=  pd:c'„v^ 


(2.9) 


with 


C'  =  KC 

771  C  771 


(2.10) 


Substituting  Equations  2.5  and  2.9  into  the  stability  criterion  W  =  with  W  =  pg(Sr- 
1)D„^  yields  a  stability  relation  in  form  similar  to  Shields  criterion  (e.g.  Raudkivi  1990) 
for  the  initiation  of  motion  of  sediment  particles 


=  (Cp-cg 


-1 


(2.11) 
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where  =  armor  specific  gravity,  g  =  acceleration  of  gravity,  and  v^.  =  critical  vertical 
velocity  at  which  armor  just  begins  to  lift.  In  terms  of  fi-om  Equation  2.2,  Equation 
2.11  becomes 


H 


(2.12) 


where  =  critical  wave  height  at  toe.  It  is  interesting  to  note  that  the  stability  number 
is  primarily  a  function  of  the  Froude  number,  v^.  /  This  formula  ties  the 

traditional  stability  relations  to  local  vertical  velocity  measurements. 

Based  on  results  of  detailed  velocity  measurements  in  the  interior  and  just 
outside  the  armor  layer,  the  vertical  velocity  gradient  was  found  to  be  proportional  to  the 
ratio  of  the  vertical  velocity  and  the  armor  diameter,  as  assumed  in  Equation  2.8.  The 
empirical  convection  coefficient  is  =  0.90  for  this  experiment.  This  is  shown  in 
Figure  2.8  for  a  group  of  experiments  summarized  in  Table  2.3.  All  experiments  listed 
in  Table  2.3  were  conducted  with  a  seaward  slope  of  1  V:2H,  D„  =  4.6  cm,  and  4  =  24 
cm.  The  velocity  values  ai  e  positive  peaks  from  the  aligned  inner  and  outer  vertical 
velocity  time  series.  In  Table  2.3,  the  velocity  gradient  Av/Ay  =  |  (v^  -  v,)/(y„  -  y,)  | ,  where 
=  outer  peak  velocity,  v,-  =  inner  peak  velocity,  y,-  =  inner  velocity  measurement 
elevation,  y^  =  outer  velocity  measurement  elevation. 


25 


Table  2.3 

Summary  of  Convection  Measurement  Experiment  Results  with  cot  6=2, 
D„  =  4.6  cm,  and  d,  =  24  cm. 


Measure¬ 

ment 

Wave 

Wave 

Location 

Height 

Period 

1 

cm 

s 

Outer 

Velocity 

Vo 

cm/s 

1/s 

Inner 

Velocity 

V, 

LvILy 

cm/s 

1/s 

near  free 
surface 


0  2  4  6  8  10  12  14 

Figure  2.8.  Vertical  variation  of  vertical  velocity  under  a  steep  wave  front  across  the 
armor  layer  as  a  function  of  wave  steepness 

The  drag  and  inertia  coefficients  can  be  more  accurately  defined  if  we 
assume  spherical  armor.  Based  on  previous  studies  of  forces  on  armor  by  Mizutani  et 
al.  (1992)  and  Torum  (1994),  reasonable  estimates  for  drag  and  inertia  coefficients  are 
Cjj  =  0.8  and  C„  =  0.4  yielding 

^  =  0-5  =  1.21)  (2.13) 

=  0.36  =  0.90)  (2.14) 

where  =  1.21  corresponds  to  a  sphere.  The  critical  vertical  velocity,  v^,  for  the 
incipient  vertical  armor  movement  reduces  to 
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(2.15) 


- ^ - =  1.2 

where  the  critical  vertical  velocity,  v^,  depends  on  the  nominal  diameter,  D„,  and  the 
specific  gravity,  only,  for  a  loose  armor  unit.  So,  if  the  vertical  velocity  exceeds  this 
critical  velocity,  motion  of  the  sphere  should  occur.  At  the  point  of  incipient  motion, 
this  critical  condition  can  be  expressed  as 

V  =  =  ^\.2D^g{S-\)  (2.16) 

Plan  4  in  Table  2.1  was  designed  to  test  the  above  criterion.  For  Plan  4,  the 
armor  layer  was  constructed  using  silicon  rubber  spheres  which  were  glued  together  and 
attached  to  an  inflexible  yet  porous  metal  mat.  The  metal  mat  was  placed  directly  on  the 
underlayer  and  fixed  to  the  flume  walls.  Several  loose  concrete  spheres  were  placed  in 
the  armor  layer  along  a  line  from  above  the  still  water  level  down  to  the  toe.  Each  two 
loose  spheres  were  separated  by  two  glued  spheres  so  that  there  was  no  interaction 
between  loose  spheres.  The  sphere  layer  of  Plan  4  was  constructed  to  have  the 
minimum  porosity  of  a  sphere  layer  of  0.33. 

For  Plan  4,  the  loose  spheres  would  not  move  under  any  conditions  unless 
they  were  slightly  raised  in  the  armor  layer.  This  was  accomplished  by  placing  a  0.5- 
cm-thick  spacer  under  each  sphere.  The  primary  effect  of  this  was  to  raise  the  porosity 
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surrounding  the  loose  sphere,  providing  a  path  for  water  motion  under  the  sphere.  The 
only  motion  observed  for  the  raised  loose  spheres  was  vertical  motion  under  the  steep 
wave  front,  following  a  slightly  elliptical  path,  and  landing  back  in  their  hole  after  the 
wave  front  passed.  For  tests  with  vertical  velocities  corresponding  to  the  critical  value, 
the  spheres  were  just  lifting  off.  For  the  larger  vertical  velocities,  the  spheres  were 
lifting  entirely  out  of  their  initial  holes,  but  settling  back  into  their  holes.  This  sphere 
motion  tmder  the  breaking  wave  is  shown  in  the  sequence  of  photographs  in  Figure  2.9. 
The  sphere  in  motion  is  just  left  of  the  black  rectangle  on  the  right  side  of  the 
photograph. 


In  the  incipient  motion  experiment,  spheres  at  a  depth  of  one-third  the  toe 
depth  were  the  most  mobile  while  spheres  at  the  stiU  water  level  were  somewhat  less 
mobile.  This  movement  corresponded  to  the  variation  of  the  vertical  velocities  in  the 
water  column  as  shown  in  Figure  2.6.  Figure  2. 10  shows  the  incipient  motion  criterion 
of  Equation  2.16  versus  wave  steepness  for  Plan  4  using  a  few  representative  points 
from  each  motion  category  as  summarized  in  Table  2.4.  For  this  figure  =  61.8  cm/s 
computed  using  Equation  2.16  with  D„  =  3.03  cm,  S,  =  2.083,  and  g  =  980.6  cm/s.  The 
dark  horizontal  line  represents  the  theoretical  incipient  motion  criteria  while  the  velocity 
measurements  are  represented  by  the  dark  dots.  Observed  movement  is  noted  for  each 
data  point.  The  vertical  gaps  between  the  lifting  group  of  points  and  the  stationary  and 
rolling  groups  occurred  because  the  vertical  velocity  increased  dramatically  under  the 
steep  breaking  wave  face.  Therefore,  it  was  difficult  to  get  a  continuous  set  of  points 
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over  the  entire  range  of  v/v^  ratios.  For  the  drag  and  inertia  coefficients  selected,  the 
incipient  motion  criteria  agrees  reasonably  well  with  the  observed  movement. 
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Figure  2.9.  Photographs  of  sphere  motion  during  a  typical  incipient  motion  experiment 
for  Plan  4 
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Figure  2.10.  Incipient  motion  criterion  (Equation  2.16  with  =  61.8  cm/s)  versus  wave 
steepness 
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Table  2.4 

Summary  of  Incipient  Motion  Experimentai  Resuits 

Sphere 

Wave 

Wave 

Wave 

Steep- 

Veiocity 

depth 

Height 

Period 

ness 

V 

V/Vj 

cm 

cm 

s 

Hr/4 

cm/s 

4.6 

1.0 

0.029 

28.3 

0.46 

13.0 

1.0 

0.083 

41.8 

0.68 

2.8 

2.0 

0.009 

11.0 

0.18 

8.4 

2.0 

0.027 

29.6 

0.48 

near  free 

3.0 

0.012 

20.2 

0.33 

surface 

3.0 

0.016 

17.2 

0.28 

d,/d,=  0.16 

14.0 

1.0 

0.090 

42.5 

0.69 

15.2 

1.5 

0.065 

67.8 

1.10 

17.0 

1.5 

0.073 

64.9 

1.05 

14.3 

2.0 

0.046 

63.8 

1.03 

13.0 

2.0 

0.042 

59.8 

0.97 

4.6 

1.0 

0.029 

8.8 

0.14 

one-third  of 

13.0 

0.083 

31.3 

0.51 

depth 

2.8 

0.009 

7.5 

0.12 

d,/d,=  0.38 

8.4 

0.027 

22.9 

0.37 

3.0 

0.012 

14.9 

0.24 

3.0 

0.016 

12.3 

0.20 
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2.7  Conclusions  from  Incipient  Motion  Study 

Experiments  on  incipient  motion  of  breakwater  armor  showed  several 
modes  of  displacement.  One  dominant  mode  was  due  to  vertical  wave  forces  which  are 
shown  to  occur  at  the  point  of  maximum  vertical  velocity  under  the  steep  wave  front.  A 
simple  relation  was  derived  assuming  a  Morison-like  wave  force  balanced  by  the  armor 
unit  submerged  weight.  The  wave  force  model  was  composed  of  drag,  due  to  the 
mflYimnm  vertical  velocities,  and  inertia,  due  to  the  vertical  convective  accelerations. 
The  maximum  vertical  convective  acceleration  is  shown  to  be  roughly  linearly  related  to 
the  square  of  the  velocity,  which  puts  the  inertial  force  term  into  the  same  form  as  the 
drag  term.  The  resulting  incipient  motion  stability  relation  is  similar  in  form  to  the 
Shields  sediment  motion  criteria.  Further,  when  expressed  as  a  traditional  stability 
number,  incipient  motion  is  shown  to  be  a  function  of  the  Froude  number,  v/{gHf\  The 
incipient  motion  criterion  shows  promise  in  predicting  the  incipient  motion  of  spheres 
for  the  conditions  tested,  but  further  experiments  are  required.  This  study  provided 
some  insight  into  how  breaking  waves  can  instigate  armor  motion  and  remove  armor 
units  from  an  intact  armor  layer.  Further  experiments  were  conducted  using  a  high- 
precision  force  transducer  to  measure  the  forces  on  exposed  spherical  armor  units. 

These  measurements  have  not  been  analyzed  completely  and  are  not  included  herein. 
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Chapter  3 


fflSTORICAL  DAMAGE  MEASUREMENT  AND  DESCRIPTION 

3.1  Damage  Modeling  Standards 

There  is  a  substantial  amount  of  literature  concerning  the  measurement  of 
damage  on  rubble-mound  coastal  structures.  Hughes  (1993)  reviewed  laboratory 
techniques  for  measuring  damage.  He  noted  three  types  of  experiments  for  accumulated 
damage:  (a)  long-duration  tests,  (b)  accumulated-storm-impacts  tests,  and  (c)  residual- 
stability  tests.  There  is  overlap  among  these  three  and  few  standards  appear  to  exist  for 
these  types  of  laboratory  studies.  Jensen  (1984)  noted  that  model  storm  duration  should 
generally  be  specified  to  provide  the  equivalent  of  8  to  10  hr  prototype.  This  is 
sufficient  if  equilibrium  damage,  where  further  waves  cause  no  additional  damage, 
occurs  in  this  time;  but  if  not,  a  subset  of  tests  should  be  conducted  to  determine  the 
ultimate  damage  level.  Hughes  stated  that  tests  should  be  repeated  at  least  two  to  four 
times  to  develop  sufficient  statistical  certainty  in  the  expected  outcome  with  more 
extensive  testing  performed  if  the  variance  is  large. 

There  are  two  dominant  methods  for  damage  measurement:  (a)  visual: 
counting  the  number  of  individual  armor  units  that  have  been  dislodged  and  moved 
more  than  one  nominal  diameter  from  their  original  location  and  (b)  profile 
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measurement:  determining  the  eroded  armor  volume  through  profiling.  A  subset  of  the 
profiling  method  of  characterizing  damage,  noted  by  Van  der  Meer  (1988)  for  dynamic 
stability,  is  through  description  of  the  profile  geometry.  Torum  et  al.  (1979)  and  Davies 
et  al.  (1994)  described  measurement  of  the  minimum  depth  of  cover,  shown  in  Figure 
3.1,  which  is  a  reduction  of  the  profile  shape  to  a  single  parameter. 


Figure  3. 1 .  Sketch  of  breakwater  profile  with  definition  of  eroded  area  and  depth  of 
cover 

3.2  Damage  Measurement  Methods 
3.2.1  Eroded  volume  method 

The  U.S.  Army  Engineer  Waterways  Experiment  Station’s  (WES’) 
historical  method  for  characterizing  damage  utilized  profiles  to  determine  the 
percentage  volume  of  stones  eroded  relative  to  the  total  volume  of  stones  in  the  active 
armor  layer.  Hudson  (1959)  used  this  volume  method.  But  the  method  apparently  is 
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not  described  in  detail  in  any  public  references.  As  such,  it  will  be  described  herein.  In 
this  volume  method,  the  active  armor  layer  was  defined  as  extending  from  the  middle  of 
the  breakwater  crest  to  one  zero-damage  wave  height  below  the  still  water  level.  The 
damage  was  determined  through  profiles  using  a  sounding  rod  with  a  circular  foot  of 
diameter  equal  to  0.56Z)„5o,  where  is  the  nominal  diameter  of  the 

median  stone  weight,  W50,  and  Ya  is  the  specific  weight  of  armor  stone.  The  sounding 
disc  size  was  determined  so  that  the  before-testing  armor  layer  thickness,  determined 
using  the  measured  profile,  coincided  with  the  theoretical  value.  The  soundings  were 
generally  obtained  on  a  horizontal  grid  spaced  evenly  at  1.5Z>„5o.  A  number  of  profiles 
along  the  breakwater  length  were  averaged  to  determine  an  average  profile.  The 
average  damaged  profile  was  subtracted  from  the  undamaged  average  profile  to  get  an 
average  eroded  area  over  the  active  region.  The  eroded  cross-sectional  area  is  defined  in 
Figure  3.1,  This  eroded  area  was  divided  by  the  total  area  of  armor  in  the  undamaged 
average  profile  to  get  a  percent  damage  D%.  Hudson’s  (1959)  zero-damage  criteria 
corresponded  to  D%  <  1  percent.  The  zero-damage  criteria  given  in  the  Shore 
Protection  Manual  (SPM  1984)  corresponds  to  D%  <  5  percent  by  the  eroded  volume 
method  or  2  percent  by  count.  The  justification  for  the  less  restrictive  zero-damage 
criteria  is  not  clear  but  evolved  over  many  years. 

The  primary  weakness  of  the  WES  eroded  volume  method  is  that,  because 
the  damage  is  only  computed  over  the  active  region,  the  damage  value  will  depend  on 


35 


the  structure  geometry.  This  method  also  provides  no  indication  of  the  severity  of 
damage  as  characterized  by  the  profile  shape. 


Jackson  (1968)  gave  damage  due  to  regular  waves  in  terms  of  the  WES 
eroded  volume  definition  of  damage.  Jackson’s  damage  values,  given  in  Table  7-9  in 
the  SPM,  have  been  widely  used  to  predict  damage  but  are  limited  to  regular  waves  with 
damage  starting  from  an  undamaged  structure.  Also,  the  damage  is  not  given  as  a 
function  of  time,  which  is  critical  for  determining  the  reliability  of  a  structure. 

Broderick  and  Ahrens  (1982)  and  Van  der  Meer  (1988)  defined  a 
dimensionless  damage  index  using  profile  data  as 


^50 


I  P- 


(3.1) 


where 

=  eroded  volume  per  unit  length  or  cross-sectional  eroded  area 
M50  =  mass  of  median  armor  unit  in  mass  distribution 

=  armor  unit  density 

where  their  variable  nomenclature  has  been  modified  to  avoid  confusion.  Broderick  and 
Ahrens  stated  that  was  calculated  from  the  profile  data  by  determining  the  difference 
between  before  and  after  testing  profiles.  The  difference  in  Broderick  and  Ahrens’ 
method  from  the  traditional  eroded  volume  method  was  that  they  nondimensionalized 
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by  the  square  of  the  nominal  stone  diameter  rather  than  the  area  in  the  undamaged  armor 
profile. 


If  the  eroded  area  over  the  entire  structure  is  used  to  compute  S  rather  than 
just  the  active  region,  the  ratio  of  damage  by  the  eroded  volume  method  to  that  of  the 
damage  index  can  be  computed  as  per  Cornett  (1995) 
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where 

-  armor  layer  thickness 
=  crest  width 

=  breakwater  crest  elevation  above  bottom 
hf  =  water  depth  at  toe 
H  =  design  wave  height 

0  =  seaside  angle  of  armor  slope  relative  to  horizontal 

Equation  3.2  assumes  the  wave  height  is  less  than  the  depth  at  the  toe. 
Cornett  noted  that  the  range  of  this  ratio  is  0.6  to  1.25  for  typical  rubble  mounds.  For 
S/D%  =  0.8,  he  notes  that  the  zero-damage  criteria  of  Z)  =  5%  corresponds  to  0  <  5  <  4. 
This  is  quite  a  broad  criterion.  The  damage  index  method  appears  to  give  a  better 
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representation  of  damage.  Broderick  and  Ahrens  and  van  der  Meer  noted  a  zero- 
damage  level  of  5  =  2. 


Thompson  and  Shuttler  (1976)  also  used  the  eroded  volume  method.  They 
used  a  profiler  with  a  foot  diameter  of  D^fjl  with  soundings  spaced  at  D„so,  which  is 
similar  to  that  used  in  the  WES  experiments.  The  structure  was  surveyed  after  each 
1,000  waves.  They  computed  the  eroded  volume  using  the  trapezoidal  rule  and  an 
average  profile.  A  damage  number  was  calculated,  assuming  a  spherical  armor 
shape,  as  the  number  of  stones  eroded  in  a  9D„^  wide  region  of  the  breakwater  section 
as  follows 


N.  = 


r^3 


nSO 


(3.3) 


which  is  equivalent  to 


54pf 
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(3.4) 


over  a  9D„5o  width  of  the  structure,  where 
=  armor  bulk  density 
Vg  =  average  eroded  volume 
=  actual  armor  density 

The  difference  between  the  method  of  Thompson  and  Shuttler  and  the 
damage  index  method  is  the  5Ap^Jnp„  in  Equation  3.4.  For  Thompson  and  Shuttler’ s 
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tests,  the  density  ratio  was  pjp^a  =1-81.  Therefore,  the  method  of  Thompson  and 
Shuttler  yielded  a  damage  index  approximately  nine  times  that  of  Broderick  and  Ahrens, 
or  the  width  of  their  structure  in  nominal  diameters.  Thompson  and  Shuttler  also 
determined  the  minimum  armor  layer  thickness  at  failure.  They  defined  failure  as  the 
point  at  which  an  area  of  exposed  underlayer  of  diameter  occurred.  Note  that  the 
minimum  armor  layer  thickness  will  not  be  zero  at  failure  because  it  is  expressed  as  a 
spatial  average  of  several  profiles. 

H.  R.  Wallingford,  Ltd.  (1990),  showed  that  Equation  3.1  yielded  very 
different  results  if  a  slightly  different  method  was  used  to  compute  the  average  eroded 
area.  The  first  method  they  used  was  that  described  for  the  WES  eroded  volume 
method,  where  an  average  profile  was  used  to  determine  an  average  eroded  area.  The 
alternative  method  was  to  sum  the  eroded  areas  from  all  profiles  in  order  to  compute  an 
average  eroded  area.  The  difference  between  the  two  methods  ranged  from  2  to  82 
percent.  In  general,  the  difference  decreased  as  the  damage  level  increased.  Note  that 
most  authors  do  not  describe  the  method  used  to  compute  damage. 

All  damage  methods  discussed  above  share  a  common  weakness,  namely 
they  compute  the  average  damage,  which  may  be  concentrated  in  one  pocket  or  spread 
out  over  several  areas.  Also,  none  of  the  methods  give  any  indication  of  the  profile 
shape,  or  more  specifically,  the  maximum  depth  of  erosion,  which  is  certainly  an 
important  parameter  for  a  multilayer  structure. 
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3.2.2  Stone  count  method 


Hedar  (1960),  Owen  and  Allsop  (1983),  Hughes  (1993),  and  Davies  et  al. 
(1994)  describe  measurement  of  damage  through  stone  counts.  Besides  visual  counting 
during  the  test,  photo  overlays  and  digital  image  processing  software  can  be  used  to 
determine  damage  by  stone  counts;  but  these  methods  are  relatively  complex  and  time 
consuming.  Stone  count  methods  are  useful  for  determining  very  low  damage  values 
but  become  inaccurate  if  more  than  a  few  stones  begin  to  move  or  if  movement  is  due  to 
sliding  rather  than  dislodgement  of  individual  stones.  Stone  count  suffers  from  the 
same  weaknesses  as  the  eroded  volume  method,  namely  that  the  spatial  concentration  of 
damage  is  generally  not  specified  and  the  maximum  depth  of  erosion  is  not  computed. 
Stone  count  is  also  somewhat  subjective. 

3.2.3  Recent  damage  measurement  methods 

Davies  et  al.  (1994)  provided  a  review  of  laboratory  techniques  for 
measuring  breakwater  profiles  and  methods  for  characterizing  damage.  They  described 
the  WES  damage  D%  (SPM  1984)  as  a  visual  displaced  stone  count  and  made  no 
reference  to  the  WES  eroded  volume  method.  They  introduced  a  technique  to  compare 
the  stone  count  method  of  damage  measurement  with  the  damage  index  method  of 
Broderick  and  Ahrens  (1982).  For  damage  measurement,  they  used  a  semiautomated 
profiler  that  measured  profiles  with  a  small  spatial  sampling  interval  by  dragging  a 
wheel  over  the  structure  face.  During  an  experiment  measuring  damage  on  a  riprap 
armor  layer,  they  computed  the  damage  index,  apparently  using  an  average  profile,  and 
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they  noted  that  stone  counts  were  more  accurate  when  only  a  few  stones  moved,  but  the 
volume  method  improved  accuracy  of  the  damage  measurement  for  advanced  damage. 
They  also  noted  that  the  depth  of  cover  was  useful  in  describing  the  degree  of 
damage.  They  noted  that,  for  an  armor  layer  thickness  of  approximately  2D„^q,  when 
the  underlayer  was  visible  through  a  hole  in  size,  and  when  d^-O, 
significant  damage  to  the  underlayer  had  occurred. 


3.3  Damage  Measurement  Experiments 

Historically,  breakwater  design  has  been  accomplished  using  an  empirical 
stability  equation,  such  as  the  Hudson  equation  (Hudson  1958,  1959)  as  shown  in 
Equation  2.1.  As  described  earlier,  for  this  equation,  Kp  is  defined  for  a  given  level  of 
performance,  typically  the  no-damage  condition  represented  by  D%  less  than  2  percent 
by  count  or  5  percent  by  volume  (SPM  1984).  This  technique  assumed  damage 
approached  an  equilibrium  level  of  D%,  where  further  regular  waves  at  the  design 
condition  induced  no  further  damage.  This  is  based  on  regular  wave  experiments  where 
damage  reaches  an  equilibrium  level  or  failure  occurs  relatively  quickly.  Hudson  (1958) 
and  Van  der  Meer  (1988)  expressed  Equation  2.1  as  a  stability  number  as  shown  in 
Equation  2.2  or 


where  A  =  5^  -1. 
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The  Hudson  equation  was  defined  for  regular  waves  but  has  been  extended 
to  irregular  waves  using  irregular  wave  physical  model  experiments.  The  SPM  suggests 
the  average  height  of  the  highest  1/10  waves,  to  replace  the  regular  wave  height  in 
Equations  2.1  and  3.5,  but  the  justification  is  unclear.  Researchers  attempted  to  relate 
regular  and  irregular  wave  effects  on  stability  in  the  1970s  with  no  definite  conclusions 
(e.g.,  Jensen  (1984)).  Recently,  Vidal  et  al.  (1991);  Vidal  et  al.  (1995);  and  Jensen  et  al. 
(1996)  emphasized  the  need  for  characterization  of  the  large  waves  in  a  random  wave 
train.  Vidal  et  al.  (1995)  noted  that  representation  of  the  irregular  time  series  by 
the  average  height  of  the  highest  100  waves,  provided  a  comparative  level  of  damage  to 
that  produced  by  regular  waves;  but  the  equivalent  statistic  of  the  form  Hy„  will  depend 
on  the  number  of  waves.  Medina  and  McDougal  (1988)  introduced  an  interesting,  albeit 
not  rigorous,  method  for  interpreting  Jackson’s  (1968)  regular  wave  damage  results 
using  a  Rayleigh  wave  height  distribution.  Their  method  incorporated  storm  duration 
into  the  equation.  In  summary,  the  regular  wave  stability  and  damage  formulations 
given  above  are  conservative  for  design;  but  a  universal  analytical  technique  for 
extending  these  relations  using  irregular  waves  and  resulting  damage  development  has 
not  been  completed. 

Thompson  and  Shuttler  (1976)  performed  both  long-term  deterioration  and 
shorter  single-storm  damage  tests  using  a  riprap-armored  embankment  with  an 
impermeable  core.  All  of  their  tests  were  restricted  to  nonbreaking  waves  in  front  of  the 
embankment,  mostly  deep  water.  The  average  run  length  was  5,000  waves,  based  on 
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mean  zero-upcrossing  wave  period  T„,  with  intermediate  surveys  every  1,000  waves. 
They  considered  5,000  waves  to  be  a  typical  storm  duration.  Several  of  their  significant 
conclusions  concerning  damage  are  paraphrased  as  follows: 


1.  The  rate  of  erosion  of  a  given  riprap  is,  as  expected, 
strongly  dependent  upon  “the  significant  wave  height” 
Hus. 

2.  The  rate  of  erosion  decreases  with  time  and  hence  the 
damage  history  curves  flatten  out.  At  the  lower  wave 
heights,  the  curves  can  become  nearly  horizontal,  giving 
an  apparently  stable  riprap  slope. 

3.  The  very  long  preliminary  tests  give  no  certainty  of  the 
riprap  eroding  to  a  totally  stable  equilibrium  state,  even 
with  low  damage  rates.  Thus,  it  is  not  safe  to  assume,  as 
is  often  done  in  regular  wave  tests,  that  a  slope  will  erode 
to  stability.  All  that  can  be  said  is  that  the  erosion  rate 
may  become  small  enough  to  be  ignored  in  practice. 

4.  The  method  of  laying  the  riprap  has  a  significant  effect 
upon  the  damage  history. 

5.  The  wave  energy  spectral  shape  as  specified  by  the 
width  parameter,  e,  does  not  affect  the  ultimate  erosion 
volume. 

6.  The  failure  criterion  requiring  a  given  area  of  exposed 
filter  layer  was  easier  to  assess  than  that  requiring  the 
observation  of  the  erosion  of  filter  material  and  gave 
erosion  volumes  at  failure  which  were  independent  of  the 
filter  grades  used. 

7.  Experimental  time  limited  the  maximum  value  of  N^, 
the  average  number  of  waves  incident  on  the  riprap,  to 
5000,  which  is  typical  of  a  storm.  In  most  cases  this  was 
too  few  waves  to  determine  whether  or  not  equilibrium 
damage  was  achieved  or  whether,  at  a  given  value  of  Hy^, 
the  slope  protection  would  eventually  fail. 
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8.  Even  in  the  limiting  case  of  waves  with  a  small 
significant  height  incident  for  a  very  long  time  on 
relatively  large  riprap,  there  will  be  a  few  rare  waves  high 
enough  to  remove  the  smallest  stones  of  the  riprap  pack 
and  hence  give  damage. 

9.  The  movement  of  the  stone  is  greatest  on  the  flatter 
slopes  although  the  net  erosion  is  small.  This  movement 
results  in  self-healing  by  the  smaller  stones. 


These  conclusions  offer  a  somewhat  different  view  of  damage  development 
than  was  accepted  at  that  time  based  on  regular  wave  experiments.  The  conclusion  that 
an  equilibrium  level  of  damage  may  not  occur  provided  motivation  to  include  damage  in 
a  stability  model. 

Using  Thompson  and  Shuttler’s  riprap  stability  data.  Van  der  Meer  (1988) 
stated  that  the  damage  rate  should  be  linear  up  to  500  to  1,000  waves  but  “for  large 
numbers  a  limit  to  the  damage  should  be  reached  (equilibrium).”  These  criteria  for  the 
relation  between  damage  and  the  number  of  waves  were  limited  to  tests  where  damage 
was  larger  than  5  =  3  after  5,000  waves  and  where  the  filter  layer  was  not  visible  after 
5,000  waves.  Van  der  Meer’s  discussion  of  damage  progression  is  limited  to  widely 
graded  armor,  which  may  not  deteriorate  in  a  manner  similar  to  uniformly  sized  stone 
armor.  This  will  be  discussed  further  in  the  next  chapter. 

3.4  Damage  Progression  Prediction 

Table  7-9  in  the  SPM  (1984),  included  herein  as  Table  3.1,  provides  a 
deterioration  model  for  armor  stability  based  on  regular  wave  data  from  model  studies 
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discussed  in  Jackson  (1968).  The  table  gives  damage  as  a  function  of  wave  height 
relative  to  the  zero-damage  wave  height  for  several  armor  types.  The  tabulated  values 
can  be  used  to  formulate  an  equation  for  armor  stone  damage.  As  an  aside,  note  that  the 
dolos  damage  progression  in  the  table  is  limited  to  stability  and  does  not  include 
breakage.  But  because  dolos  movement  causes  breakage,  this  damage  progression  may 
not  be  conservative. 


Table  3.1 

Jackson  (1968)  Breakwater  Damage  from 
Table  7.9  of  SPM  (1984) 

Percent 

Rough 

Smooth 

Damage 

Stone 

Stone 

D% 

HIH,, 

2  to  5 

1.00 

1.00 

5  to  10 

1.08 

1.08 

lOtolS 

1.19 

1.14 

15  to  20 

1.27 

1.20 

20  to  30 

1.37 

1.29 

30  to  40 

1.47 

1.41 

40  to  50 

1.56 

1.54 

A  simple  empirical  model  for  the  best  fit  line  through  the  data  points  in  the 
table  for  rough  stone  damage  is  given  by 


D%  =/ 


H 


H 


do 


2.79  exp 


V  "doj 


-  18.08 


(3.6) 
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The  experiments  supporting  Table  3.1  were  conducted  with  regular  waves,  so  measured 
damage  approached  an  equilibrium  after  a  relatively  short  duration.  Therefore  damage 
is  not  a  function  of  storm  duration  in  Equation  3.6.  Figure  3.2  shows  data  from  Table 
3.1  and  damage  given  by  Equation  3.6.  The  relation  utilizes  H,  the  regular  wave  height, 
and  Hjg,  the  no-damage  regular  wave  height  corresponding  to  0  to  5  percent  damage  by 
the  eroded  volume  method.  This  damage  relation  shows  how  damage  varies  with  wave 
height  but  it  assumes  starting  with  an  undamaged  structure  and  damage  approaching  an 
equilibrium  level. 
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Figure  3.2.  Damage,  characterized  as  eroded  volume,  as  a  function  of  monochromatic 
wave  height  for  angular  stone  (SPM  1984) 


Van  der  Meer  (1988)  reanalyzed  data  from  Thompson  and  Shuttler  (1976) 
and  conducted  a  number  of  additional  experiments  and  found  the  damage  index  S  to  be 
related  to  the  number  of  waves  V,,.,  the  significant  wave  height  and  the  mean  wave 
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period  T„.  Van  der  Meer  noted  that  S  =  2  in  Equation  3.1  provides  a  good  estimate  of 
the  initiation  of  damage  and  that  failure  occurred  when  5  =  8  for  structure  slopes  of 
IV:  1.5H  and  1V:2.0H,  where  failure  was  defined  as  exposure  of  the  underlayer  through 
a  hole  D„5o  in  diameter.  Van  der  Meer  used  test  durations  of  1,000  and  3,000  average 
wave  periods.  Virtually  all  of  his  tests  were  performed  with  nonbreaking  waves.  Van 
der  Meer  performed  eight  tests  with  depth  limited  waves  and  eight  more  where  perhaps 
only  the  highest  waves  in  the  distribution  were  depth  limited.  Based  on  these  data,  he 
proposed  breaking- wave-induced  stability  after  V,,,  waves  could  be  determined  by  the 
following  equations. 


For  plunging  waves; 


H. 


2% 


AD 


nSO 


=  1.4(6.2)P®‘* 


p-0.5 

9m 


(3.7) 


and  for  surging  waves: 


H. 


2% 


AD. 


nSO 


1_4P-0.13[ 


0.2  _ 

V/C^C 


(3.8) 


where 

^2%  =  wave  height  exceeded  by  2  percent  of  waves  in  wave  height 

probability  distribution,  where  Van  der  Meer  assumed  jF/2%  =\AH^ 
to  obtain  Equations  3.7  and  3.8  from  his  formulations  for 
nonbreaking  waves  at  toe  of  structure 

A  =  (5^-1)  with  Sg  =  specific  gravity  of  armor  unit 

P  =  permeability  coefficient 

0  =  structure  slope  angle 
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=  tan0  /  (H/Lj^  =  surf  similarity  or  Iribarren  parameter 
H,  =  significant  wave  height  which  is  equal  to  if, ,3 
4  =  gT^Jiln) 

=  mean  wave  period 

Van  der  Meer’s  breaking  wave  tests  were  performed  with  a  traditional 
multi-layer  rubble  mound  where  the  notional  permeability  was  P  =  0.5.  This 
permeability  is  an  empirical  parameter  without  any  regard  to  the  flow  throughout  porous 
media.  Using  this  value  of  P,  Equations  3.7  and  3.8  suggest  that  S  is  approximately 
proportional  to  and  Because  the  number  of  waves  is  defined  as  N^=  t/T„, 
where  is  the  total  run  time,  these  formulations  indicate  that  damage  increases  with  the 
square  root  of  time.  But  because  the  Iribarren  parameter  is  raised  to  a  negative  constant 
in  Equation  3.7  but  raised  to  a  positive  power  of  P  in  Equation  3.8,  the  effect  of  wave 
period  is  not  clear.  It  is  clear  that  damage  progression  is  very  sensitive  to  wave  height. 

Equations  3.7  and  3.8  are  of  limited  practical  use  for  depth-limited  breaking 
waves  primarily  because  the  supporting  breaking  wave  experiment  was  extremely 
limited  in  scope.  These  equations  are  essentially  the  same  as  Van  der  Meer’s 
nonbreaking  wave  equations  except  the  Rayleigh  wave  height  distribution  assumption 
of  =  H2%I\A  has  been  substituted.  Van  der  Meer  included  the  structure  slope  in  the 
Iribarren  or  surf  similarity  parameter  (Battjes  1974),  but  the  beach  slope  is  critical  for 
depth-limited  breaking  waves.  Van  der  Meer  conducted  experiments  for  only  one  beach 
slope  Md  used  a  veiy  narrow  range  of  wave  periods.  Therefore,  the  effect  of  varying 
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beach  slope  or  breaker  type  was  not  determined.  For  breaking  waves,  as  the  beach  slope 
steepens,  the  wave  breaking  becomes  more  vigorous  (Battjes  1974)  and  the  slamming 
forces  increase  dramatically  as  the  waves  start  plunging  to  collapsing  (Bruun  1985). 

Van  der  Meer’s  spilling  breakers  were  not  very  severe  with  respect  to  stability  relative 
to  plunging  or  collapsing  breakers. 

Equations  3.7  and  3.8  do  not  indicate  decreasing  stability  as  the  depth  to 
wave  length  ratio  decreases  and  the  wave  breaking  becomes  more  severe.  Carver  and 
Wright  (1991)  showed  that  stability  decreases  dramatically  with  decreasing  relative 
depth.  The  minimum  stability  condition  occurs  where  the  wave  breaking  is 
characterized  by  plunging  to  collapsing  breakers  at  the  toe  of  the  strucmre.  Figure  3.3 
shows  some  of  Carver  and  Wright’s  data  replotted  with  Hudson  stability  coefficient, 
in  Equation  2.1,  as  a  function  of  relative  depth.  Here  the  water  depth  h,  at  the  structure 
toe  is  normalized  by  the  local  wave  length  Lp  based  on  the  peak  spectral  period 
computed  using  linear  wave  theory.  It  is  clear  that  the  Hudson  stability  coefficient 
decreases  with  decreasing  relative  depth.  This  is  a  reflection  of  the  severity  and  location 
of  the  breaking  wave  and  the  resulting  wave  forces.  As  was  shown  in  Chapter  2,  the 
critical  wave  forcing  for  incipient  motion  was  uplift  occurring  under  the  steep  wave  face 
for  collapsing  to  plunging  breakers.  For  depth-limited  conditions,  as  the  water  depth 
decreases  relative  to  the  wave  length,  the  wave  face  steepens,  and  vertical  forces  at  the 
steep  wave  face  are  able  to  loosen  and  mobilize  the  stones.  Therefore,  for  depth-limited 


49 


waves,  the  critical  condition  for  stability  occurs  for  relatively  long  waves  that  plunge  on 
the  structure. 

A  minor  difficulty  with  Equation  3.7  and  3.8  is  a  potential  source  of 
confusion  because  the  stability  number  is  given  in  terms  of  H2%,  but  the  Mbarren 
parameter  is  in  terms  of  Finally,  all  of  the  previous  damage  empirical  equations 
listed  above  are  limited  to  constant  incident  wave  conditions  and  water  level  as  well  as 
the  initial  condition  of  an  undamaged  structure.  These  equations  are  intended  to  predict 
damage  for  a  single  design  storm.  Therefore,  they  cannot  be  used  to  predict  damage 
development  over  the  life  of  a  structure  as  is  required  in  a  life-cycle  cost  analysis. 

3,5  Variability  in  Stability  Results 

Carver  and  Wright  (1991)  conducted  irregular  wave  stability  experiments 
primarily  intended  to  determine  random  variations  of  damage  initiation  due  to  varying 
armor  placement  and  due  to  variations  in  wave  time  series  realizations  with  constant 
spectral  parameters.  Their  tests  progressed  to  low  damage  levels,  varying  up  to  7.7 
percent  displacement,  by  count,  of  armor.  The  wave  conditions  were  assumed  to  be 
constant  for  a  given  test.  But  they  did  vary  the  wave  height  and  period  between 
structure  rebuilds  if  the  damage  had  not  progressed  far  enough.  This  testing  strategy  is 
not  typical  and  their  interpretation  of  the  data  assumed  damage  caused  by  one  wave 
condition  would  be  independent  of  the  damage  caused  by  the  following  wave  condition. 
They  made  note  of  the  uncertainty  in  stability  of  traditional  stone-armored  breakwaters. 
In  Figure  3.3,  for  each  value  of  relative  depth,  the  stability  coefficient  is  shown  to 
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scatter  widely.  They  concluded  that  the  uncertainty  was  due  to  differing  construction 
techniques  from  test  to  test  and  different  random  number  seeds  being  used  to  generate 
the  spectra. 

Font  (1968)  noted  that  the  armor  placement  technique  is  important  for 
determining  when  initiation  of  damage  occurs,  but  is  less  relevant  for  advanced  damage, 
which  is  similar  to  conclusions  stated  above  by  Thompson  and  Shuttler.  It  is  anticipated 
that  the  scatter  in  damage  will  be  large  for  advanced  damage  primarily  because  of  the 
random  nature  of  the  interaction  between  irregularly  shaped  stones.  Melby  and  Mlakar 
(1997)  showed  how  Carver  and  Wright’s  data  could  be  used  to  compute  the  reliability 
of  a  breakwater  with  respect  to  the  zero-damage  condition  of  armor  stability  and  how 
Van  der  Meer’s  relation  could  be  used  to  compute  reliability  of  a  damaged  mound.  But 
little  work  has  been  done  on  computing  the  reliability  of  a  damaged  rubble  mound 
through  simulated  damage  using  historical  storm  data. 
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Figure  3.3.  Hudson  stability  coefficient  versus  relative  depth  for  angular  armor  stone 
exposed  to  irregular  waves  (Carver  and  Wright  1991) 
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Chapter  4 


DETERIORATION  EXPERIMENT 


4.1  Overview 

The  deficiencies  of  previous  damage  experiments  were  noted  in  the 
preceding  chapters.  Because  accurate  quantification  of  damage  is  required  for  a  risk 
analysis,  a  small-scale  laboratoiy  experiment  was  conducted  to  quantify  breakwater 
deterioration.  The  initial  experiment  utilized  a  traditional  trapezoidal  multilayer  rubble- 
mound  armored  with  uniform-sized  stone  which  was  exposed  to  depth-limited  irregular 
waves.  The  objectives  of  the  experiment  were  as  follows: 

a.  Determine  deterioration  rates  for  the  structure  under  long-term  exposure  to 
storm  conditions. 

b.  Develop  and  evaluate  methods  for  damage  characterization. 

c.  Queintify  development  of  damage  with  changing  storm  parameters,  e.g.,  water 
level,  wave  height. 

d.  Quantify  profile  development  as  related  to  damage. 

e.  Quantify  alongshore  variability  in  damage. 

f.  Evaluate  and  standardize  methods  for  breakwater  profiling. 
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4.2  Experimental  Setup 
4.2.1  Flume  and  structure 

The  experiment  was  accomplished  using  two  identical  side-by-side 
structures  constructed  in  a  flume  measuring  61.1m  long,  1.52  m  wide,  and  2.0  m  high 
(Figure  4.1).  The  structure  cross  sections  are  shown  in  Figures  4.2  and  4.3.  Each  0.76- 
m-wide  structure  was  constructed  on  a  flat  bottom.  The  offshore  concrete  beach  slope 
was  1  V:20H.  This  steepness  was  selected  to  realistically  produce  severely  breaking 
waves  in  the  vicinity  of  the  structure  toe. 


Wave  Generator 

Offshore 


Nearshore  Wave 
Gages ' 


Test 

Structures  Wave 

Absorber 


Figure  4.1  Flume  profile  for  damage  experiments 


Figure  4.2  Model  structure  cross  section  for  damage  experiments 


54 


Figure  4.3.  Photograph  of  model  structures  before  testing 

In  armor  stability  experiments,  breakwater  crest  heights  are  usually  set  high 
enough  to  prevent  overtopping  from  influencing  the  stability,  unless  overtopping  is  the 
focus  of  the  study  or  the  breakwater  is  low  crested.  But  the  composite  downslope 
weight  of  upslope  armor  units  generally  contributes  to  the  stability  of  the  armor  layer 
and  increases  as  the  slope  lengthens;  although  under  some  conditions,  longer  armor 
layers  can  have  a  tendency  to  buckle  or  slide.  Just  considering  stone  stability,  for  a 
given  wave  height,  taller  structures  are  typically  more  stable,  with  respect  to  main  armor 
stability,  than  short  structures.  But,  in  the  prototype,  the  cost  of  the  breakwater  is  kept 
to  a  minimum  by  setting  the  breakwater  crest  height  such  that  minor  overtopping  is 
allowed  for  a  design  wave  condition.  Thus,  experimenting  with  non-overtopped 
structures  that  have  high  crest  heights  typically  will  not  be  conservative  with  respect  to 
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damage  development.  For  the  experiment  discussed  herein,  the  structure  crest  height 
above  the  toe  of  30.5  cm  corresponded  to  an  elevation  of  approximately  one  significant 
wave  height  above  the  still  water  level  for  the  larger  wave  conditions  tested.  This  is 
typical  of  prototype  structure  crest  heights.  This  crest  height  was  specified  so  that 
results  from  this  experiment  could  be  used  for  analysis  of  prototype  breakwaters  with 
higher  crest  heights  because  the  damage  measured  with  the  shorter  structure  should  be 
relatively  greater  than  that  of  a  taller  structure.  The  effect  of  wave  overtopping  on 
stability  will  be  discussed  in  more  detail  later  in  this  chapter. 

4.2.2  Experimental  test  conditions 

The  test  conditions  for  the  initial  study  are  summarized  in  Table  4, 1 ,  Note 
that  wave  period  variation  and  armor  gradation  variation  experiments  are  not  listed  in 
this  table,  but  were  added,  and  are  discussed  in  Chapter  7.  The  waves  had  relatively 
long  periods  in  this  initial  experiment,  as  this  produces  the  lowest  stability  (Figure  3.3). 
The  depth-to-wavelength  ratio  h,/Lp  at  the  stmcture  toe  was  approximately  0.07  for  this 
initial  experiment,  where  was  computed  using  linear  wave  theory  and  the  spectral 
peak  period.  This  shallow-water  breaking  wave  is  also  typical  of  design  conditions  on 
United  States  coastlines. 

4.2.3  Armor  layer  and  underlayer 

The  mass  distributions  for  armor  and  underlayers  for  the  initial  test  series 
are  shown  in  Figure  4.4.  The  armor  layer  consisted  of  uniformly  hand-sized  stones  with 
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Table  4.1 

Experimental  Conditions  for  Initial  Damage  Experiment 

Parameter 

Range 

Water  depth  at  toe,  h, 

11.9-15.8  cm 

Wave  height  at  toe 

Variable 

Peak  wave  period, 

2.48  -  2.59  sec 

Seaward  beach  slope,  cot  a 

20 

Structure  slope,  cot  6 

2 

Structure  crest  height  above  bottom, 

30.5  cm 

Nominal  armor  stone  diameter, 

3.64  cm 

Median  armor  stone  mass, 

128  g 

Stone  density,  p„ 

2.66  g/cm^ 

Mjo  =  128  g,  while  the  underlayer  was  sieve  sized  passing  1.59-cm  and  retained  on 
1.27-cm  sieves.  The  ratio  of  armor  median  mass  to  underlayer  median  mass  was 
(M5o)a  /  (M^oX  =  25,  which  was  large  relative  to  that  typically  used  in  the  prototype.  The 
gradation  of  the  armor  layer  for  this  initial  series  was  very  narrow  with  =  1.05, 

where  Dgj  is  the  nominal  diameter  corresponding  to  the  85-percent  exceedance 
probability  in  stone  mass  and  is  the  nominal  stone  diameter  corresponding  to  the  15- 
percent  exceedance  probability.  The  gradation  of  the  underlayer  was  wider  with 
Dgj/Dis  =  1.44.  The  core  material  was  sieve  sized,  passing  No.  3  (0.67-cm)  and  retained 
on  No.  4  (0.47-cm)  sieves.  While  not  impermeable,  this  core  material  was  sized  small 
enough  to  prevent  strong  flow  within  the  core.  Equations  3.7  and  3.8  indicate  that  lower 
porosity  in  the  underlayer  and  core  results  in  lower  stability.  Therefore,  the  underlayer 
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and  core  were  selected  to  be  at  the  lower  end  of  the  range  of  material  that  is  used  in  the 
prototype. 


Figure  4.4.  Armor  and  underlayer  stone  mass  distributions  for  initial  experiment 

The  amount  of  material  in  the  armor  layer  is  determined  using  the  packing 
density  equation,  given  in  the  SPM  as 


A 

where 

=  number  of  armor  units 
A  =  unit  area  of  breakwater  slope  to  be  armored 
n  =  number  of  armor  layer  thicknesses 
=  layer  coefficient 
=  armor  layer  porosity  in  percentage 
W5o=  median  armor  weight 


(4.1) 
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It  will  be  convenient  to  refer  to  damage  in  terms  of  the  nominal  stone 
diameter,  so  the  packing  density  equation  can  be  rearranged  to  be  compatible  as 


A  7)2 


(4.2) 


where  <{)  is  the  packing  density  coefficient.  The  stone  armor  layer  was  placed  in  a 
traditional  two-layer  thickness  at  an  average  packing  density  coefficient  of  (j)  =  1.2, 
corresponding  to  n  =  2,  =  1.0,  and  40  percent  as  recommended  in  the  SPM. 

The  core  was  placed  by  shovel,  troweled  to  grade,  and  washed  in  place  to 
naturally  pack  tight.  This  process  simulates  the  natural  washing  action  of  waves  during 
the  construction  sequence.  Similarly,  the  underlayer  was  dumped  onto  the  sloping  core 
and  lightly  troweled  to  grade.  Each  armor  stone  was  placed  individually  by  hand, 
simulating  crane  placement  as  closely  as  possible.  The  armor  stone  placement  during 
the  experiment  followed  WES  guidelines  for  random  hand  placement  in  the  laboratory 
as  follows: 

a.  Stones  were  lowered  vertically  into  position. 

b.  Stones  were  placed  on  the  slope  rather  than  dropped. 

c.  Stones  were  placed  so  that  they  touched  their  neighbors. 

d.  Stones  could  not  be  pushed  into  a  hole. 

e.  Particular  orientation  of  stones  or  placement  patterns  were  avoided. 
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4.2.4  Model-Prototype  similitude 

Many  authors  have  discussed  similitude  requirements  for  breakwater 
stability  including  Dai  and  Kamel  (1969),  Hudson  et  al.  (1979),  Sakakiyama  and  Kajima 
(1992),  and  Hughes  (1993).  As  such,  only  a  brief  review  will  be  provided  herein.  In 
order  for  damage  measurements  in  the  physical  model  to  be  scalable,  similitude  must  be 
maintained,  where  the  ratios  of  the  dominant  physical  forces  in  the  model  are  the  same 
as  they  are  in  the  prototype.  Because  there  is  no  mathematical  relationship  that 
describes  armor  damage,  dimensional  and  inspectional  analysis  and  the  resulting 
empirical  relations  are  relied  on  to  establish  similitude  requirements  (Hudson  et  al. 

1979,  Hughes  1993). 

The  first  requirement  for  similitude  is  that  the  model  be  geometrically  and 
kinematically  undistorted.  Therefore,  the  size  and  shape  of  the  breakwater,  armor  and 
underlayer  stone,  surface  characteristics  of  the  stone,  and  the  height  and  length  of  the 
waves  must  all  be  undistorted.  As  described  in  Chapter  2,  the  wave  forces  on  the  armor 
units  act  to  cause  damage  while  the  armor  unit  self-weight  and  inter-unit  friction  act  to 
prevent  armor  movement.  Also  shown  in  Chapter  2,  wave  forces  are  typically 
decomposed  into  fluid  drag,  fluid  inertia,  and  buoyancy.  The  two  important  force  ratios 
that  include  these  forces  are  the  Froude  number  and  the  Reynolds  number  (e.g.  Hudson 
et  al.  1979).  The  Froude  number  is  the  square  root  of  the  ratio  of  the  inertial  and 
gravitational  forces  or  Fr  =  uligtj^  where  u  is  the  fluid  velocity,  g  is  the  gravitational 
acceleration,  and  /  is  a  characteristic  dimension  of  an  armor  unit  or  wave  height.  Froude 
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similitude  demands  that  Fr„  =  Fr^,  where  m  and  p  refer  to  model  and  prototype, 
respectively.  The  Reynolds  number  is  the  ratio  of  the  inertial  forces  to  viscous  forces 
Re  =  uUv,  where  v  is  the  fluid  kinematic  viscosity.  So,  in  addition  to  geometric 
similitude,  complete  similarity  in  damage  requires  Froude  and  Reynolds  similitude  or 


Fr. 


p  _ 


Fr_ 


=  1, 


Re. 


(4.3) 


In  stability  models,  Froude  and  geometric  similimde  are  typically  maintained  leading  to 
the  relationship  between  length  and  time  scales  of 


N,  -  W.  -  ^  (4.4) 

where  Nj,  and  Ni^  are  the  time,  velocity,  and  length  scales,  respectively.  Under 
Froude  similitude,  in  order  to  achieve  Reynolds  similitude,  the  following  relationship 
must  be  satisfied. 


JV.  ■  nP  (4.5) 

where  is  the  viscosity  scale  ratio.  Therefore,  it  is  not  practically  possible  to  attain 
both  Reynolds  and  Froude  similitude  with  the  use  of  water  in  small-scale  stability 
experiments.  Water  is  always  used  in  stability  models,  yielding  some  Reynolds  scale 
effects.  Thus  the  viscous  forces  are  relatively  stronger  in  the  model  than  in  the 
prototype. 
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In  order  to  maintain  reasonably  close  similarity  between  model  and 
prototype,  armor  stability  models  are  usually  constructed  at  a  scale  large  enough  to 
avoid  significant  Reynolds  scale  effects  by  ensuring  that  the  flow  around  the  armor 
remains  turbulent.  Dai  and  Kamel  (1969)  suggested  that  Re  =  /v  >  3  x  10“* 

will  prevent  Reynolds  scale  effects.  Their  results  were  based  on  large-  and  small-scale 
flume  tests  using  regular  waves.  Van  der  Meer  (1988)  stated  that  the  lower  limit  of 
Reynolds  numbers  should  be  in  the  range 

Jig  =  >  1x10^-4x10''  (4.6) 

v 

to  prevent  Reynolds  scale  effects.  This  guidance  is  presently  accepted  practice.  In  the 
experiment  discussed  herein,  the  Reynolds  number  range  is  3.1  x  ICC  <  Re  <  4.0  x  10'' 
which  satisfies  the  requirements  suggested  by  both  Dai  and  Kamel  and  van  der  Meer. 

4.3  Wave  Generation  and  Measurement  in  Initial  Experiment 
4.3.1  Wave  generation 

Waves  were  generated  based  on  the  Texel,  Marsen,  and  Arsloe  (TMA) 
spectrum  (Bouws  et  al.  1985).  The  deterministic  spectral  amplitude  and  random  phase 
method  was  used  to  synthesize  the  time  series  that  was  used  to  drive  the  piston-type 
wave  board.  Long  (1986)  describes  the  computer  program  used  to  synthesize  the  time 
series.  The  software  creates  a  variable  spectral  bandwidth  Ao)„  by  the  method  of  Goda 
(1970),  which  allows  closer  spectral  lines  near  the  peak  frequency.  For  this  method,  a 
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sequence  of  random  phases  0„  are  specified  over  the  interval  from  0  to  27t.  Then  the 
Fourier  amplitudes  are  found  from  the  relation 


=  V2S^(a>„)Ao)„  (4.7) 

given  a  discrete  target  spectrum  as  a  function  of  angular  frequency  («)„.  The  final  step 
is  to  determine  the  discrete  time  series  by  inverse  Fourier  transforming  the  amplitudes 
and  phases.  This  method  provides  a  very  close  match  with  the  target  peak  frequency 
and  the  spectral  shape. 

For  these  initial  test  series,  two  wave  board  drive  signals  were  generated 
corresponding  to  two  water  depths.  The  design  conditions  for  wave  generation  are 
listed  in  Table  4.2.  Note  that  the  gain  is  generator  specific  but  provides  a  crude  judge  of 
the  energy  levels  because  gain  is  roughly  linearly  related  to  wave  height. 


Table  4.2 

Wave  Generation  in  initial  Damage 

1  Experiment 

Wave 

Generator 

Depth 

cm 

Peak 

Period 

sec 

Gain 

as  percent  of 
maximum 

1 

112.7 

2.50 

50 

2 

112.7 

2.50 

70 

3 

112.7 

2.50 

90 

4 

116.5 

2.50 

50 

5 

116.5 

2.50 

70 

6 

116.5 

2.50 

90 
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The  random  number  seed  value  for  the  random  phase  generation  was 
different  for  the  two  depths.  These  two  signals  were  scaled  to  achieve  three  wave 
heights  for  each  water  depth  by  adjusting  the  gain  on  the  output  amplifiers  for  the  wave 
board  command  signal  as  indicated  in  Table  4.2.  The  resulting  six  time  series  of  the 
command  signals  sent  to  the  wave  board  are  shown  in  Figure  A.l  in  Appendix  A.  Note 
that  short  segments  at  the  beginning  and  end  of  each  series  have  been  cut  off  in  the 
analysis.  The  spectra  for  these  six  signals  are  shown  in  Figures  A.2  and  A.3,  for  the  two 
water  depths.  The  wave  board  was  commanded  at  a  rate  of  20  Hz,  and  ramps  of  5  sec 
each  were  placed  at  the  beginning  and  end  of  the  command  time  series  to  gradually  start 
and  stop  the  wave  generator. 

4.3.2  Wave  measurement 

Waves  were  measured  using  capacitance-type  wave  gages.  Two  arrays  of 
three  gages,  offshore  and  nearshore,  were  used  as  shown  in  Figure  4.1.  Nearshore  gage 
layout  is  shown  in  Figure  4.5.  Considering  the  largest  wave  measured  herein,  this  gage 
location  approximately  corresponds  to  the  recommendation  by  Goda  (1985)  where  the 
design  breaking  wave  is  determined  seaward  of  the  toe.  This  is  the  travel  distance 
of  large  breaking  waves.  For  the  toe  depth  oih,=  11.9  cm,  the  depth  at  the  most 
nearshore  gage  was  16.5  cm,  and  the  depth  at  the  offshore  array  was  1 12.7  cm.  For  h,  = 
15.8  cm,  the  depth  at  the  most  nearshore  gage  was  20.4  cm  and  the  depth  at  the  offshore 
array  was  1 16.5  cm. 
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Landward  Middle  Seaward 

Gage  Gage  Gage 


Figure  4.5.  Hume  profile  near  structure  with  nearshore  wave  gage  layout 

The  gages  were  calibrated  using  the  following  automated  technique.  Each 
gage  was  attached  to  a  Jordan  servo-controller  and  motor.  The  wave  gage  data 
acquisition,  control,  and  signal  conditioning  system  consisted  of  a  VAX  PDP-1 1 
microcomputer  connected  to  12-bit  differential  input  analog-to-digital  converters, 
digital-to-analog  converters,  and  analog  filters.  A  program  on  the  PDP-1 1  sent  signals 
to  the  Jordan  controllers  sending  each  through  10  calibration  stops  twice,  covering  a 
range  of  motion  in  excess  of  the  largest  wave  expected  in  the  irregular  time  series.  At 
each  stop,  the  voltage  from  each  gage  was  recorded.  For  each  gage,  a  calibration  curve 
was  automatically  tested  for  fit  through  the  calibration  points  using  linear,  quadratic,  or 
cubic  spline  regression.  The  calibration  program  automatically  chose  the  simplest 
technique  with  less  than  0.25  mm  deviation  from  the  regression  curve,  which 
corresponds  to  approximately  0.1 -percent  maximum  deviation  over  the  full  calibration 
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range.  This  was  the  maximum  error  allowed  and  the  error  varied  because  the  calibration 
range  was  different  for  each  gage. 

Measured  nearshore  water  surface  oscillation  time  series  are  shown  in 
Figures  A.4  through  A.6  for  waves  of  three  significant  heights  at  /i,  =  1 1.9  cm.  Figures 
A.7  through  A.9  show  wave  gage  time  series  for  waves  of  three  significant  heights  at  h, 
=  15.8  cm.  Incident  and  reflected  waves  were  resolved  using  these  time  series  from  the 
three-gage  arrays  using  the  technique  of  Goda  and  Suzuki  (1976).  The  technique  was 
modified  as  per  Kobayshi,  Cox,  and  Wuranto  (1990)  to  determine  the  incident  and 
reflected  wave  spectra  and  time  series  at  the  shallowest  gage.  Time  series  parameters 
were  determined  using  the  zero-upcrossing  method.  At  frequencies  containing  little 
wave  energy,  the  reflection  analysis  technique  can  yield  poor  estimates  due  to  low 
signal-to-noise  ratios.  These  regions  are  indicated  by  low  coherency  in  the  cross 
correlation  between  gage  pairs.  For  this  study,  the  cutoff  frequencies  were  established 
to  maintain  the  coherence  above  0.3  in  the  cross  correlation  between  gage  pairs.  The 
analysis  technique  was  not  sensitive  to  coherence  cutoffs  higher  than  0.3. 

The  resolved  incident  and  reflected  spectra  for  time  series  from  the 
nearshore  and  offshore  arrays  are  shown  in  Figures  A.IO  through  A.12  for  =  1 1.9  cm 
and  in  Figures  A.13  through  A.15  for  h,  =  15.8  cm.  Figures  A.IO  and  A.13  show  an 
increase  in  peak  energy  density  due  to  wave  shoaling  between  the  offshore  and 
nearshore  gages  and  an  increase  in  wave  energy  for  low  (/<  0.35  Hz)  and  high  (f  >  l.O 
Hz)  frequencies.  These  figures  show  that  wave  energy  is  being  transferred  due  to 
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nonlinear  wave  interactions  from  the  frequency  band  0.35  -1.0  Hz  to  frequency  bands 
/<  0.35  Hz  and/  >1.0  Hz.  Figures  A.1 1,  A.  12,  A.  14,  and  A.  15,  for  the  larger  wave 
heights,  show  a  similar  transfer  of  energy  but  also  show  a  decrease  in  peak  wave  energy 
due  to  wave  breaking.  It  is  clear  from  these  figures  that  the  wave  energy  decrease  due  to 
wave  breaking  increased  dramatically  as  the  generated  wave  energy  increased.  It  is  also 
interesting  that  the  measured  reflected  wave  energy  was  much  higher  nearshore  than 
offshore.  This  trend  is  primarily  because  wave  breaking  and  bottom  friction  reduce  the 
incident  wave  energy  as  waves  progress  shoreward  while  the  reflected  waves  undergo 
little  energy  change  as  they  progress  offshore.  This  is  consistent  with  the  measured 
cross-shore  variations  of  wave  reflection  from  beaches  (Baquerizo  et  al.  1997). 

Tables  4.3  and  4.4  list  nearshore  incident-wave  characteristics  for  the  six 
wave  trains,  each  of  15-min  duration,  measured  with  the  structure  in  place.  Table  4.3 
lists  the  spectral  analysis  parameters,  and  Table  4.4  lists  the  time  series  parameters,  both 
at  the  location  of  the  shallowest  gage  that  was  0.91  m  seaward  of  the  toe.  The  values 
were  computed  using  the  analysis  method  discussed  above  with  the  three  nearshore 
gages.  In  Table  4.3,  the  mean  period  is  defined  by  the  relation 

,  J/S-K/W 

-  (4.8) 

\Uf)df 

0 

where  S,fj)  is  the  energy  spectral  density  of  the  incident  wave  and /is  the  frequency. 

The  spectral  significant  wave  height  is  defined  here  as  with  m  =  zero 


67 


moment  of  the  incident  wave  spectram.  The  average  reflection  coefficient  is  defined  as 
R  =  [(mj/m^y^  with  (mj^  =  zero  moment  of  the  reflected  wave  spectrum.  The  change 
of  wave  reflection  with  damage  measurement  was  negligible  in  this  experiment.  This  is 
consistent  with  the  reflection  measurements  made  by  Smith  et  al.  (1992)  on  a  horizontal 
bottom  in  deeper  water  in  front  of  a  1  V:3H  slope.  The  local  wave  length  L„  and  the 
deep  water  wave  length  =  gTjIln  were  both  computed  using  linear  wave  theory  and 

the  mean  period.  The  Iribarren  parameter  in  Table  4.3  is  defined  as 

=  tan  and  in  Table  4.4  as  =  tan  where  tan  a  =  1/20  is 

the  nearshore  beach  slope  and  in  Table  4.4. 


Table  4.3 

Nearshore  Wave  Statistics  for  incident  Spectra  and  Related 
Parameters  with  Structure  in  Place 

■ 

Refiec- 

Relative 

Relative 

Spectral 

Iribarren 

■ 

Toe 

Wave 

Mean 

Peak 

tion 

Toe 

Toe 

Para- 

Depth 

Height 

Period 

Period 

Coeff. 

Depth 

Depth 

meter 

■ 

b, 

cm 

^mo 

cm 

Tm 

sec 

T, 

sec 

R 

h,/L„ 

bf  7^oni 

1 

11.9 

9.78 

1.72 

2.48 

0.59 

0.066 

0.026 

0.34 

2 

11.9 

12.4 

1.70 

2.48 

0.60 

0.067 

0.027 

0.30 

3 

11.9 

14.2 

1.72 

2.48 

0.60 

0.066 

0.026 

029 

D 

15.8 

10.5 

1.72 

2.59 

0.61 

0.076 

0.034 

0.33 

15.8 

13.6 

1.69 

2.59 

0.62 

0.077 

0.035 

0.29 

15.8 

15.8 

1.68 

2.59 

0.61 

0.078 

0.036 

0.26 
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Table  4.4 

Nearshore  Wave  Statistics  for  incident  Time  Series  and  Related 

1  Parameters  with  Structure  in  Place 

Wave 

Toe 

Depth 

ft. 

cm 

Wave 

fteight 

^1/3 

cm 

Period 

Tn, 

sec 

stability 

Number 

Iribarren 

Parameter 

^tom 

Period 

sec 

Wave 

Height 

cm 

Period 

^lAIO 

sec 

Wave 

Height 

WlAtO 

cm 

Stability 

Coefficient 

fti^f(Hino) 

1 

11.9 

9.38 

1.76 

1.6 

0.36 

2.16 

12.8 

2.28 

11.5 

3.4 

2 

11.9 

11.6 

1.69 

1.9 

0.31 

2.13 

15.4 

2.40 

13.8 

6.0 

3 

11.9 

13.2 

1.74 

2.2 

0.30 

2.22 

17.3 

2.53 

15.7 

8.8 

D 

15.8 

10.1 

1.73 

1.7 

0.34 

2.15 

14.3 

2.31 

12.7 

4.7 

5 

15.8 

13.0 

1.67 

2.2 

0.29 

2.16 

17.8 

2.39 

16.0 

9.3 

6 

15.8 

14.9 

1.66 

2.5 

0.27 

2.20 

19.3 

2.47 

18.0 

13.2 

For  a  Rayleigh  distribution  of  wave  heights,  =  1.27  and  = 

1.40.  Computing  these  ratios  using  the  values  from  Table  4.4,  -  1.19-1.26  and 

=  1.30-1.42,  indicating  the  slight  reduction  of  wave  heights  due  to  wave 
breaking. 

Traditionally,  the  incident  waves  are  measured  without  the  structure  in  place 
(Hughes  1993).  For  this  study,  the  waves  in  the  flume  were  measured  both  with  and 
without  the  structure  in  place.  For  the  condition  with  no  structure  in  place,  an  array  of 
three  gages  was  placed  with  its  centroid  at  the  location  of  the  structure  toe  for 
calibration  purposes.  The  wave  gage  array  was  moved  to  the  location  shown  in  Figure 
4.5  when  the  structure  was  in  place.  Wave  heights  measured  without  the  structure  in 
place  were  smaller  than  those  measured  with  the  structure  in  place,  and  the  difference  in 
ranged  from  6  to  24  percent.  The  difference  between  the  two  sets  of  measurements 
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was  due  to  the  fact  that  the  wave  generation  system  did  not  allow  absorption  of  reflected 
and  re-reflected  waves  within  the  flume;  so  wave  energy  built  up  in  the  flume  during 
testing.  It  should  be  noted  that  systems  that  absorb  reflected  energy  are  only  partially 
effective,  and  so  it  is  virtually  impossible  to  prevent  wave  energy  increase  due  to  re¬ 
reflected  waves  in  stability  tests. 

The  wave  energy  in  this  experiment  reached  an  equilibrium  state  quickly 
during  each  ran  due  to  the  relatively  long  wave  period.  A  simple  analysis  was  done  to 
determine  the  rate  of  increase  of  wave  energy.  For  this  analysis,  the  wave  time  series 
were  divided  into  eight  segments,  each  of 4,096  points  or  3.41  min,  except  the  first 
segment  which  was  2,048  points.  Both  the  first  and  second  segments  started  at  the 
400th  data  point  or  t  =  20  s  and  the  final  segment  ended  at  the  16,784th  data  point  or  r  = 
14.0  min.  The  segments  overlapped  by  2,048  data  points,  and  each  segment  represented 
approximately  120  mean  wave  periods,  except  the  first  segment,  which  was  about  60 
wave  periods  long.  For  each  segment,  the  incident  was  calculated  using  the 
technique  described  above.  Figures  4.6  and  4.7  show  the  wave  height  variation  for  h,  = 

1 1.9  cm  and  h,  =  15.8  cm,  respectively.  In  these  figures,  the  ratio  of  the  incident 
segment  wave  height  to  the  overall  wave  height  is  shown  as  a  function  of 

normalized  duration  where  is  the  duration  to  the  center  of  the  segment  and  = 
15  min  is  the  total  ran  duration.  As  can  be  seen  from  the  figures,  the  wave  height 
reached  an  equilibrium  very  early  during  each  ran.  Based  on  these  observations,  the 
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wave  heights  measured  with  the  structure  in  place  were  used,  as  they  were  more 
accurate  than  those  measured  without  the  structure  in  place. 


Figure  4.6.  Wave  height  variation  throughout  Waves  1, 2,  and  3  at  1 1.9  cm  toe  depth 


The  wave  reflection  analysis  as  used  assumed  a  flat  sea  bottom  and  that 
linear  wave  theory  was  valid.  But  the  bottom  of  the  flume  sloped  at  1V:20H  and  the 
waves  were  clearly  nonlinear.  An  analysis  was  done  to  determine  the  sensitivity  of  the 
wave  analysis  method  to  selection  of  gage  location  and  therefore  beach  slope.  Using  the 
depth  at  the  center  of  the  gage  array  in  the  reflection  analysis,  rather  than  the  depth  at 
Gage  1,  resulted  in  no  difference  for  2-percent  difference  for  and  7-percent 
difference  for  The  difference  in  depths  between  Gage  1  and  the  array  center  was 
12  percent.  and  were  more  stable  than  in  this  analysis.  H,  and  for  Gage 
1  closest  to  the  toe  were  approximately  the  same  for  repeated  tests.  This  is  important,  as 
the  damage  was  shown  to  be  approximately  proportional  to  the  fifth  power  of  the  wave 
height  in  the  previous  chapter,  exaggerating  any  errors  in  wave  height  measurement. 
Therefore,  the  more  stable  values  of  and  are  relied  on  in  the  analysis  that 
follows.  Comparison  of  measured  wave  heights  with  breaker  index  curves,  such  as 
Figure  7-2  of  the  SPM,  showed  that  the  maximum  Hj/jo  would  just  break  at  the  mean 
period  in  the  depth  of  the  shallowest  wave  gage,  approximately  1  m  from  the  structure 
toe. 


The  Iribarren  parameter  or  surf  similarity  parameter  (Battjes  1974)  is 

defined  as 


I  = 


tan0 


(4.9) 
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where  =  gT^/lTi  and  tan  0  is  either  the  structure  slope  or  the  beach  slope,  depending 
on  where  the  wave  breaks.  The  wave  height  H  is  usually  either  or  and  the 
period  is  the  mean  period  for  armor  stability  (Van  der  Meer  1988).  The  Iribarren 
parameter  can  be  used  to  characterize  the  type  of  wave  breaking  with  ^  <  0.3 
corresponding  to  spilling  waves,  0.3  ^  2  to  plunging  waves,  2  ^  3  to  collapsing 

waves,  and  ^  s  3  to  surging  waves.  Van  der  Meer  (1988)  and  others  have  used  the 
structure  slope  in  Equation  4.9.  This  is  because  they  were  primarily  testing  with 
nonbreaking  waves  on  the  beach.  But  for  the  experiment  described  herein,  the  wave 
heights  were  depth  limited.  So  the  beach  slope  was  the  critical  parameter  dictating  the 
type  of  wave  breaking  near  the  structure  and  was  therefore  used  to  compute  It  is  clear 

from  Tables  4.3  and  4.4  that,  for  these  tests,  the  most  severe  waves  were  primarily 
plunging.  This  corresponds  to  the  design  condition  for  most  United  States  shoreline 
applications.  The  relative  depths  show  that  the  wave  condition  selected  is  near  the 
worst  case  for  stability  according  to  Figure  3.3. 

The  stability  numbers  shown  in  Table  4.4  were  in  the  range  =  1. 6-2.5. 

This  is  well  within  the  range  of  =  1-4  for  conventional  breakwaters  suggested  by  Van 
der  Meer  (1988).  On  the  other  hand,  stability  coefficients  in  Table  4.4  based  on  //,/jo  are 
in  the  range  3.4-13.2  in  comparison  to  Kp  =  2.0  recommended  in  the  SPM  for  breaking 
waves  on  rough  angular  stone  armor  placed  on  a  trunk.  This  indicates  that  the 
breakwaters  in  this  experiment  will  incur  some  damage. 
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4.4  Effect  of  Overtopping  on  Stability 

The  effect  of  wave  overtopping  on  armor  stability  was  investigated  by  Van 
der  Meer  and  Daemen  (1994).  They  gave  a  reduction  factor  for  stable  armor  size  as 


= - 1 -  for  0<Rp<  0.052 

1.25 -4.8i?/ 


(4.10) 


with 


'  h.N 


op 

2n 


IT" 

sT, 


(4.11) 


where  is  the  freeboard  or  structure  crest  height  above  the  still  water  level.  There  is  an 
odd  mix  of  wave  height  statistics  in  Equation  4.11.  Substituting  the  conditions  in  the 
experiment  discussed  herein  yields  R*  of  0.063  -  0.076  for  the  toe  depth  of  hj=  1 1.9  cm 
and  0.047  -  0.058  for  h,  =  15.8  cm.  Therefore,  the  1 1.9-cm  depth  is  outside  the  range  of 
applicability  of  Equations  4.10  and  4.1 1.  Actually,  wave  overtopping  was  negligible  in 
the  shallow  depth  of  /z,  =  1 1 .9  cm.  For  h,- 15. S  cm,  the  range  of  stable  armor  size 
reduction  coefficients  is  0.98  <k„<  1.0.  Therefore,  wave  overtopping  had  negligible 
effect  on  stability  of  the  armor  layer. 


4.5  Armor  Profiler  and  Profile  Measuring  Technique 

A  structure  profiler  similar  in  concept  to  that  used  by  Davies  et  al.  (1994) 
was  constructed.  The  profiler  design  and  construction  are  described  in  detail  in 
Winkelman  (1998),  who  participated  in  the  planning  and  setup  of  the  project. 

The  profiler,  shown  in  Figures  4.8  and  4.9,  consisted  of  eight  aluminum  profiling  arms. 
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laterally  spaced  at  5  cm,  which  pivoted  about  a  point  on  a  sliding  carriage.  A  rotating 
spheroid  with  diameter  approximately  D„so  at  the  end  of  each  profiling  arm  followed  the 
structure  surface  as  the  carriage  passed  over  the  structure,  from  just  landward  of  the 
crest  to  just  past  the  seaward  toe. 


Breakwater  Toe 
Sea  Side 


Figure  4.8.  Plan  (top)  and  profile  (bottom)  views  of  breakwater  profiler 
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DETERIORATION 

TESTS 


TESTING 


Figure  4.9.  Photograph  of  breakwater  profiler 
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Note  that  acoustic  and  laser  profiling  instruments  were  available  as 
alternatives  to  the  mechanical  system  used,  but  the  signal  returns  off  the  irregular  stone 
surfaces  were  poor  resulting  in  unacceptably  high  data  analysis  requirements  and  low 
accuracy.  Also,  these  backscatter-type  sensors  were  generally  an  order  of  magnitude 
more  costly.  The  spheroid  diameter  was  chosen  to  coincide  with  that  of  Davies  et  al. 
(1994)  so  that  results  could  be  compared  and  was  nearly  twice  the  diameter  of  the 
sounding  disc  in  the  WES  eroded  volume  method.  In  the  present  case,  the  sphere 
provides  varying  resolution  of  voids  depending  on  the  depth  of  penetration  between 
stones  and  so  corresponds  to  a  range  of  sounding  disc  sizes.  No  attempt  is  made  herein 
to  relate  the  sphere  size  to  the  WES  sounding  disc  size  since  there  exist  no  damage 
progression  data  from  the  historical  records  that  are  directly  comparable  with  data 
obtained  for  the  present  study.  The  profiler  arm  pivot  point  was  at  a  fixed  elevation  just 
above  the  crest  of  the  structure  with  the  profiler  arms  extending  back  toward  the 
structure.  The  angle  of  an  arm  at  any  point  during  the  profiling  process  was  determined 
through  recorded  voltages  from  rotational  potentiometers  attached  to  the  pivot  on  the 
carriage.  The  location  of  the  carriage  was  determined  using  a  linear  potentiometer  fixed 
to  the  flume  wall. 

It  was  required  that  the  profiler  arms  be  just  heavy  enough  to  continuously 
follow  the  structure  without  plowing,  or  moving  stones.  Another  requirement  was  that 
the  structure  be  profiled  without  lowering  the  water  level  between  runs.  Therefore,  the 
hinge  point,  rod  length,  rod  and  ball  weight,  and  counter  weight  were  chosen  to  provide 
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an  optimal  balance  over  the  length  of  the  slope,  where  the  profiling  rods  were  out  of 
water  at  the  crest  and  partially  submerged  at  the  toe. 

The  power  input  to,  and  analog  output  from,  each  of  the  potentiometers  was 
controlled  by  a  personal  computer  (PC)-based  data  acquisition  system.  Analog-to- 
digital  data  conversion  and  signal  conditioning  and  recording  were  performed  using  an 
Optim  MEGADAC  5000.  The  MEGADAC  had  4  MB  of  internal  memory  and  could 
sample  up  to  250  kHz  from  20  channels.  The  heart  of  the  MEGADAC  consisted  of 
AD682SH  modules  that  included  12-bit  differential  input  analog-to-digital  converters 
and  programmable  Bessel  filters.  The  MEGADAC  also  included  power  amplifiers  to 
provide  power  to  the  potentiometers.  The  digital  voltages  were  recorded  to  a  hard  disk 
within  the  PC  at  20  Hz.  The  speed  of  the  carriage  over  the  structure  varied,  so  this 
sampling  rate  produced  a  variable  spatial  sampling  resolution.  But,  in  general,  the 
horizontal  cross-shore  sampling  resolution  was  approximately  one  sample  per 
millimeter.  Optim,  Inc.  TCS  software  was  utilized  to  control  data  acquisition. 

Prior  to  each  day’s  tests,  the  profiler  was  calibrated  by  first  recording  the 
rotational  potentiometer  voltages  with  the  rods  positioned  horizontally  on  a  jig  that  was 
set  using  high-resolution  tiltometers.  Then  the  rods  were  swung  through  a  known  angle 
and  the  voltage  differences  recorded.  The  resulting  voltage-to-angle  ratios  provided  a 
daily  specific  calibration  coefficient  for  each  potentiometer.  Prior  to  each  profile,  the 
output  voltage  for  each  potentiometer  was  set  to  zero  on  the  fixed  bar  just  above  the 
breakwater  crest. 
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The  elevations  of  each  sphere  bottom  above  the  toe,  Za,  and  horizontal 
distances  from  an  arbitraiy  starting  point,  x^,  were  determined  using  the  potentiometer 
measurements.  The  profiler  sphere  horizontal  coordinate  was  computed  as  + 
L(l-cos  (J)),  relative  to  =  0  with  (|)  =  0,  and  the  profiler  sphere  vertical  elevation  above 
bottom  was  computed  as  Za  =  Zp-  Lsin  <})  -  Dll,  where  L  =  length  of  the  profiler  arm 
from  pivot  point  to  center  of  sphere,  (})  =  angle  of  the  arm  from  horizontal,  D  =  sphere 
diameter,  =  transverse  distance  the  carriage  moved,  and  Zp  =  elevation  of  the  arm 
pivot  point.  Figure  4.8  shows  a  definition  sketch  for  this  transformation,  where  the 
initial  carriage  position  =  0  corresponds  to  the  initial  sphere  position  =  0  with  (|)  = 
0.  The  measured  sphere  elevations  were  averaged  over  a  0.5-cm  interval  in  the  cross¬ 
shore  direction  to  remove  small  variations  in  elevations  resulting  from  minor  stone 
settlements  that  occurred  during  the  tests.  This  is  an  essential  step  as  this  “noise”  will 
cause  a  bias  in  the  eroded  area  calculation. 

For  each  arm,  the  eroded  area  was  calculated  using  a  Simpson’s  Rule 
numerical  integration  of  the  damaged  profile  below  the  undamaged  profile.  The  eroded 
depth  was  determined  as  the  maximum  slope-normal  distance  between  undamaged  and 
damaged  profiles.  This  slope-normal  distance  was  calculated  using  a  numerical 
technique  beginning  at  each  point  on  the  eroded  profile  and  stepping  along  a  line  normal 
to  the  slope  of  cot  d  =  2  until  the  undamaged  profile  was  crossed.  Then  the  point  at 
which  the  crossing  occurred  was  interpolated  to  obtain  the  slope-normal  distance.  The 
slope-normal  minimum  depth  of  cover  at  a  given  time  t  was  calculated  as  the  minimnni 
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vertical  difference  between  eroded  profile  and  underlayer  multiplied  by  the  cosine  of  the 
structure  slope  angle  or 

=  min[z/j:,t) -z„(-*,t=O)]cos0  (4.12) 

where  min  indicates  the  minimum  value  with  respect  to  x,  Za{x,t)  is  the  damaged  profile 
elevation  at  location  ;c  and  time  t,  and  z„(x,/=0)  is  the  underlayer  elevation  at  ;c  before 
testing  began.  The  numerical  interpolation  technique  used  to  compute  the  eroded  depth 
was  not  required  for  the  depth  of  cover  as  the  underlayer  was  smooth  relative  to  the 
profiler  probe  size  in  comparison  to  the  undamaged  armor  profile. 
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Chapter  5 


DAMAGE  MEASUREMENTS 


5.1  Overview 

In  all,  five  test  series  were  conducted  with  two  structures  in  the  flume  within 
this  initial  experiment.  The  first  test  was  not  repeated.  The  other  four  series  included 
two  originals  and  two  repeats;  so  three  unique  series  were  conducted  and  will  be 
described  in  Chapters  5  and  6:  Series  A',  B\  and  C'.  As  noted  earlier,  additional  test 
series  were  added  to  investigate  wave  period  and  armor  gradation  and  they  are  discussed 
in  Chapter  7.  Table  5.1  describes  the  test  sequences  for  Series  A',  B’,  and  C.  Tables  B1 
through  B3  in  Appendix  B  provide  a  more  detailed  summary  of  the  Series  A',  B',  and  C. 


Table  5.1 

Summary  of  Test  Series  for  Initial  Damage  Experiment 

Test 

Series 

Test 

Type 

Order  of  Wave 
Conditions 
(Tables  2  and  3) 

Water-Level 

Order 

Test  Duration, 
hr 

B 

Deterioration  to 
failure 

1,2,  3,  4,  5,6 

Low  -  High 

28.5 

B' 

1,2, 3,  5,6 

Low  -  High 

8.5 

C 

Storm  ordering 

4,  5,  6,  2,  3 

High  -  Low 

9.0 
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For  Series  A',  lasting  28.5  hr,  the  general  experimental  strategy  was 
designed  to  provide  an  indication  of  the  long-term  deterioration  of  a  structure.  The 
structure  was  exposed  to  waves  until  failure,  where  the  underlayer  was  visible  through  a 
hole  in  the  armor  layer  at  least  in  diameter.  For  this  series,  the  experiment  began  at 
the  lower  water  depth  and  progressed  through  each  of  the  three  wave  heights  listed  in 
Tables  4.3  and  4.4,  reaching  an  apparent  equilibrium  (visually  observed)  at  each.  Then 
the  water  level  was  raised,  and  each  of  the  three  wave  conditions  at  the  greater  depth 
was  run  until  an  apparent  equilibrium  was  achieved.  The  structures  were  exposed  to  a 
total  of  approximately  60,000  waves  for  Series  A'.  This  is  equivalent  to  100  to  250  hr  of 
accumulated  prototype  storms  for  mean  wave  periods  ranging  from  7  to  15  sec.  The 
number  of  waves  was  computed  by  dividing  the  run  length  by  the  mean  period  of  the 
incident  wave  time  series. 

For  Series  B'  and  C',  the  testing  strategy  was  intended  to  simulate  damage 
from  a  sequence  of  individual  storms.  So  for  B'  and  C',  each  wave- water  level  condition 
was  run  for  approximately  4,200  waves,  representing  roughly  10-hr  prototype  storms. 

For  B’  and  C,  each  sequence  was  run  twice.  And  because  there  were  two  structures  in 
the  flume,  this  testing  provided  an  original  and  three  repeats  to  determine  statistical 
variabilities.  The  total  run  length  for  B'  and  C  was  8.5  and  9  hr  or  approximately 
18,000  waves.  Damage  did  not  progress  to  failure  for  Series  B'  and  C'.  In  Series  B',  the 
three  wave  heights  at  the  low  water  level  were  run  first  (Waves  1, 2,  3),  followed  by  the 
two  highest  wave  heights  at  the  high  water  level  (Waves  5, 6).  Series  C  was  similar 
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except  the  high-water  condition  was  run  first  followed  by  the  low-water  condition 
(Waves  4, 5, 6, 2, 3).  Throughout  the  experiment,  the  waves  were  run  in  15  min  bursts, 
with  the  water  completely  settling  between  wave  bursts.  Profiles  were  measured  after 
each  two  of  these  bursts,  every  30  min  of  waves. 

5.2  Damage  and  Eroded  Profile  Parameters 

In  Chapter  3,  damage  was  characterized  by  the  nondimensional  eroded  area 
5  defined  by  Equation  3.1,  and  this  parameter  is  used  to  describe  damage  herein.  For  a 
damage  analysis,  the  shape  of  the  eroded  portion  of  the  slope  may  be  characterized  by 
three  parameters  of  primary  engineering  interest  shown  in  Figure  5.1:  cover  depth 
eroded  depth  d^,  and  eroded  length  4-  The  cover  depth,  the  minimum  remaining  depth 
of  the  cover  layer  along  a  profile,  was  shown  by  Torum  et  al.  (1979)  and  Davies  et  al. 
(1994)  to  be  a  useful  parameter  to  characterize  the  reserve  capacity  of  the  armor  layer. 
For  this  experiment,  the  nondimensional  cover  depth  was  computed  as  the  minimum 
depth  of  armor  layer  remaining  for  each  measured  profile  after  each  30  min  of  waves 
and  normalized  by  the  nominal  stone  diameter  as 


where  was  computed  using  Equation  4.12  and  C  is  the  normalized  minimum  cover 
depth. 
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Armor  Layer 


Figure  5. 1 .  Sketch  of  breakwater  profile  with  definition  of  damage  parameters 

In  this  study,  eroded  depth  and  eroded  length  are  introduced  to  describe  the 
shape  of  the  eroded  area.  They  are  parameters  of  primary  engineering  interest, 
providing  a  description  of  extent  of  damage  normal  to  and  along  the  slope,  respectively. 
The  eroded  depth  is  defined  as  the  maximum  difference  between  undamaged  and 
damaged  profiles  measured  normal  to  the  slope  and  is  normalized  as 


(5.2) 


The  eroded  depth  is  useful  because  it  gives  an  indication  of  progress  toward  failure, 
particularly  if  failure  is  defined  by  exposure  of  the  underlayer.  The  sum  of  eroded  depth 
and  cover  depth  is  approximately  equal  to  the  initial  cover  layer  thickness  4 

=  (£  +  Q  (5.3) 
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which  is  only  approximate  because  the  maximum  and  the  minimum  may  not  occur 
at  the  same  location  due  to  variability  in  the  original  layer  thickness. 


The  eroded  length  ^  is  defined  herein  as 


de 


(5.4) 


and  the  normalized  eroded  length  is  given  by 


L  = 


_  25 


(5.5) 


Here  the  dimensional  eroded  length  4  is  not  measured  directly,  but  is  a  derived  value.  It 
is  roughly  equivalent  to  the  surface  length  of  the  eroded  region  along  a  profile.  As  seen 
in  Figure  5.1,  the  eroded  length  characterizes  the  along-profile  length  of  the  roughly 
triangular  eroded  area  and  describes  the  along-slope  extent  of  damage. 


5.3  Probability  Density  Functions  of  Damage,  Eroded  Depth,  and  Cover  Depth 

The  number  of  damage  profiles  measured  at  a  given  time  during  damage 
progression  is  denoted  by  Np  in  the  following,  where  Ap  =  16  for  Series  A'  and  Np  =  32 
for  Series  B’  and  C'.  The  Np  values  for  5,  E,  and  C  are  analyzed  statistically  to  reduce 
the  number  of  measured  values  and  separate  the  variability  along  the  structure  from  the 
temporal  variation  associated  with  damage  progression.  First,  the  mean  and  standard 
deviation  of  the  Np  values  of  S,  E,  and  C  are  calculated.  The  standard  deviations  o^ 
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and  Oc,  of  S,  E,  and  C  indicate  the  variability  of  S,  E,  and  C  among  the  profiles.  The 
means  and  standard  deviations  of  S,  E,  and  C  vary  with  damage  progression. 

The  Np  values  of  S,  E,  and  C  are  normalized  as 


S*  = 


5-5 


£■*  = 


E  -  E 


C  -  C 


(5.6) 


where  the  means,  indicated  by  overbars,  and  the  standard  deviations,  o,  of  S ,  E ,  and 
are  zero  and  unity,  respectively.  As  an  example,  the  mean  and  the  standard  deviation 
for  Series  A'  after  each  30  min  of  waves  were  computed  as 


Np 

Np  „-.i 


16 


(5.7) 


= 


Np- 16 


(5.8) 


Nptf, 


where  S„  is  the  individual  measured  value  of  5.  The  probability  density  functions  P(S ), 
P(E*),  and  P(C*)  are  calculated  using  the  Np  values  of  5*,  £*,  and  C‘,  respectively,  where 
the  range  between  the  maximum  and  minimum  values  is  divided  into  eight  bins  of 
constant  spacing.  Admittedly,  eight  bins  do  not  yield  a  fine  resolution,  but  the  number 
of  data  points  for  this  estimation  is  limited  to  =  16  or  32. 


In  summary,  the  Np  damaged  profiles  measured  every  0.5  hr  yield  one  data 
set  consisting  of  means  and  standard  deviations  of  5,  E,  and  C  and  the  probability 
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density  functions  of  S*,  E*,  and  C*.  The  number  of  data  sets  is  57  for  Series  A'  lasting 
28.5  hr,  17  for  Series  B'  lasting  8.5  hr,  and  18  for  Series  C'  lasting  9.0  hr.  For  the 
undamaged  profile  at  the  beginning  of  each  test  series,  S  =  0  and  E  =  0,  which  yield 
*^0  =  0,  =  0,  .^  =  0,  and  =  0,  where  the  subscript  zero  indicates  initial  values.  For 

an  undamaged  structure,  the  probability  density  functions  of  5  and  E  are  the  delta 
functions  with  a  spike  at  S  =  0  and  E  =  0,  but  S*  and  E*  defined  in  Equation  5.6  cannot 
be  computed  because  a^o  =  0  and  =  0.  On  the  other  hand,  the  normalized  cover 
depth  C  for  the  undamaged  profile  is  a  statistical  variable  related  to  the  placement  of  the 
armor  layer.  The  Np  undamaged  profiles  are  analyzed  in  the  same  way  to  obtain  the 
mean  Q,  the  standard  deviation  o^o,  and  the  probability  density  function  P(C*)  at  the 
beginning  of  each  test  series. 

The  probability  density  functions  of  5*,  £*,  and  C*  exhibit  large  scatters  but 
do  not  indicate  any  specific  variations  with  damage  progression  and  among  the  three 
test  series.  All  92  data  sets  of  P(5*)  and  P{E*)  are  plotted  together  in  Figures  5.2  and 
5.3,  respectively,  whereas  all  95  data  sets  of  P{C*)  are  plotted  in  Figure  5.4.  The 
standard  normal  distribution  given  by 

P(y)  =  exp(-0.5y==)/v/27r  (5.9) 

with  y  =  S*,  E*,  or  C*  is  also  shown  in  each  figure  for  reference.  Note  that  a  fit  to  the 
distribution  was  not  attempted  because  each  data  set  consisted  of  too  few  points  to 
determine  the  fit  error.  The  data  points  in  these  figures  scatter  about  the  normal 
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distribution.  The  approximate  ranges  of  the  data  points  in  these  figures  may  be  chosen 
as 


-2.7  <  5*<  3  ;  -2.7  <  E*<  2.7  ;  -2.7  <  C*<  2.8  (5.10) 

which  are  in  the  range  from  -3  to  3.  The  coirespondirig  ranges  of  5,  E,  ^d  C  can  be 
found  using  Equations  5.6  and  5.10.  This  approximate  procedure  allows  one  to  estimate 
the  lower  and  upper  limits  of  S,  E,  and  C  in  terms  of  their  means  and  standard 
deviations,  which  vary  with  damage  progression. 


S‘ 


•  Series  A'  =  Series  B’  »  Series  C - Normal 


Figure  5.2.  Probability  density  function  for  normalized  damage  S*  for  Series  A',  B',  and 
C 
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•  Series  A*  =  Series  B'  »  Series  C’ - Normal 


Figure  5.3.  Probability  density  function  for  normalized  eroded  depths*  for  Series .4' 
B\  and  C' 


C* 


Series  A'  °  Series  B'  »  Series  C - Normal 


Figure  5.4. 
andC' 


Probability  density  function  for  normalized  cover  depth  C*  for  Series  A B\ 
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5.4  Temporal  Damage  Development 

The  means  and  standard  deviations  of  5,  E,  and  C  are  used  hereafter  to 
examine  damage  progression.  Figure  5.5  shows  the  structures  prior  to  Series  A'. 

Figures  5.6  through  5.8  show  the  measured  profiles  at  the  beginning,  midway  through, 
and  at  the  end  of  Series  A'.  In  Figure  5.6,  one  profile  and  the  average  of  all  eight 
profiles  on  one  of  the  structures  at  the  beginning  of  Series  A'  are  shown.  In  Figure  5.7, 
the  average  profile  is  shown  along  with  a  photo  of  the  structures  at  the  midway  point  of 
the  series.  This  figure  corresponds  to  an  average  damage  level  oi  S=  6.5.  It  may  be 
seen  from  the  photo  that  the  top  stones  in  the  cover  layer  have  been  eroded  sporadically 
over  a  wide  region  near  the  still  water  level.  Figure  5.8  shows  one  profile  and  the 
average  profile  at  the  end  of  Series  A',  where  the  underlayer  is  exposed  through  a  hole  of 
Z)„5o  diameter.  Figure  5.9  shows  photos  of  the  structures  following  Series  A'.  It  is  clear 
that  most  of  the  armor  layer  has  been  mobilized. 
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Figure  5.5.  Photographs  of  undamaged  structures  prior  to  Series  A 
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EtEVflTIOH  FROM  lOH  (cn)  ElEUflllOM  FROM  TOE  ten) 


Cage  1  underlayer,  danaged,  and  undanaged  profiles 

Profile  Data  Fron  File: 

aprl2_18.0a 

-  Underlayer  Profile 


Average  underlayer,  danaged,  and  undanaged  profiles 

Profile  Data  Fron  File: 
aprl2_18.B8 


Underlayer  Profile 


Figure  5.6.  Profiles  at  beginning  of  Series  A',  damage  level  of  5*  =  0 
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ELEUATION  FROH  TOE  (cm) 


Average  underlayer,  danaged.  and  undanaged  profiles 

Profile  Data  From  File: 
apr01_10.04 


-  Underlayer  Profile 


0  Z0  40  60  BB  lee  120 

HORIZONTAL  SCALE  (cn) 


Figure  5.7.  Profiles  midway  through  Series  A',  damage  level  of  S=  6.5 


93 


ELEVATION  FRON  TOE  (cn)  ELEVATION  FRON  TOE  (cn) 


Gage  1  underlayer,  danaged,  and  undanaged  profiles 

Profile  Data  Fron  File: 
aprGl_31.13 


Underlayer  Profile 


Average  underlayer,  danaged.  and  undanaged  profiles 

Profile  Data  Fron  File: 
aprei_31.13 


Underlayer  Profile 


Figure  5.8.  Profiles  following  completion  of  Series  A',  dainage  level  of  5'=  12.8 
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Figure  5.9.  Photographs  of  structures  following  completion  of  Series  A',  damage  level 
ofS=  12.8 
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Figure  5.10  shows  the  number  of  waves  iV,,  versus  the  measured  average 
damage  5 for  Series  A'  for  A/p  =  16.  Here,  the  cumulative  number  of  waves  is 
obtained  using  the  mean  wave  period  T„  from  Table  4.4  and  duration  of  Waves  1-6 
listed  in  Table  B 1  in  Appendix  B.  Also  shown  are  the  alongshore  variability  of  damage 
as  one  standard  deviation  above  and  below  S.  For  Series  A',  waves  were  run  at  the 
initial  depth  until  the  structure  appeared  to  stabilize,  which  occurred  after  approximately 
22,000  waves.  At  that  point,  the  water  depth  was  raised.  As  noted  in  Table  4.4,  the 
stability  number  ranged  from  1.6  to  2.2  for  the  shallower  depth  and  from  1.7  to  2.5  for 
the  deeper  depth.  In  this  figure,  the  significant  wave  heights  H,  and  associated  stability 
numbers  N,  are  listed  across  the  top  of  the  figure,  while  the  water  depths  at  the  toe  are 
listed  across  the  bottom  of  the  figure. 


Figure  5.10.  Number  of  waves  versus  mean  damage  ±  one  standard  deviation  for  Series 
A' 
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Figure  5.10  shows  that  the  rate  of  damage  increased  with  wave  height,  and 
damage  became  more  scattered  as  the  damage  increased.  Although  the  wave  heights 
were  reduced  when  the  depth  was  increased  approximately  midway  through  the  run,  the 
damage  continued  to  increase.  This  was  due  to  the  fact  that,  when  the  depth  was 
increased,  several  stones  higher  on  the  structure  (at  the  new  still  water  level)  were 
displaced  down  into  the  damage  area,  which  was  centered  at  the  lower  still  water  level. 
This  can  be  seen  in  Figure  5.10  as  small  jumps  in  eroded  area  after  the  depth  change. 
These  small  increments  of  damage  approximately  correspond  to  one  stone  displaced. 

Although  the  structure  visually  appeared  to  stabilize  at  various  times  during 
Series  A',  Figure  5.10  indicates  that  the  structure  never  really  stabilized  or  reached 
equilibrium,  where  further  waves  produced  no  additional  damage.  Thompson  and 
Shuttler  (1976)  noted  this  for  their  riprap  tests.  It  is  also  clear  that  an  extraordinarily 
large  number  of  waves  was  required  to  induce  failure.  For  Series  A',  approximately 
60,000  waves  at  significant  wave  heights  appreciably  greater  than  the  no-damage  wave 
height  were  required  to  fail  the  structure.  Again,  failure  was  defined  as  exposure  of  the 
underlayer  through  a  hole  at  least  in  diameter.  A  higher  mean  damage  value  at 
failure  was  recorded  for  Series  A'  than  was  noted  by  Van  der  Meer  (1988).  Van  der 
Meer  quoted  a  failure  damage  value  of  5=  8,  whereas  5=  13  at  failure  for  Series  A'.  It 
was  noted  during  the  series  that  the  waves  had  a  difficult  time  moving  the  second  layer 
of  stones.  This  led  to  high  damage  values  at  failure,  with  most  of  the  first  layer  of 
stones  moved  before  the  second  layer  began  moving.  Most  other  recent  damage  testing 
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has  been  done  with  riprap,  which  consists  of  stone  having  a  wide  mass  distribution 
(Thompson  and  Shuttler  1976;  Davies  et  al.  1994).  Riprap  may  not  exhibit  this 
secondary  armoring  because  the  small  stones  in  the  mass  distribution  are  continuously 
moving.  Deterioration  of  a  riprap  layer  is  examined  in  Chapter  7. 

The  sawtooth  nature  of  the  damage  curves  was  caused  by  partial  self- 
healing.  If  a  stone  moved  out  of  its  location  and  started  downslope,  the  next  profile 
would  show  a  hole  and  relatively  increased  damage.  But  during  the  following  wave 
burst,  the  hole  would  fill  in  through  adjustment  of  surrounding  stone  so  that  the  next 
profile  would  show  a  partially  healed  structure.  This  healing  process  can  be  clearly  seen 
in  Figure  5.10  just  after  the  water  level  was  raised.  Over  time,  the  structure  became 
more  and  more  loosened,  to  the  point  that  it  could  no  longer  heal  itself  because  the 
surrounding  stones  were  also  displaced. 

Series  B'  and  C  were  similar  to  Series  A'  in  that  the  wave  and  water  level 
sequences  were  similar;  but  the  durations  and  ordering  of  the  storms  were  different  for 
B'andC'.  The  results  of  Series  S' and  C' are  tabulated  in  Appendix  B.  Mean  damage 
versus  number  of  waves  for  the  storm  sequences  of  Series  B’  and  C  are  shown  in 
Figures  5.1 1  and  5.12.  Also  shown  are  the  variability  as  a  standard  deviation  above  and 
below  the  mean.  For  these  series,  the  32  profiles  were  used  to  compute  each  data  point, 
corresponding  to  the  four  repeats  with  eight  profiles  each.  The  characteristics  of  armor 
deterioration  for  Series  B'  were  similar  to  Series  A'  because  the  order  of  the  storms  was 
the  same  for  both  series.  But  for  Series  C',  Waves  4-6  at  the  high-water  level  increased 
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damage  initially.  Waves  2  and  3  at  the  lower  water  followed  but  caused  little  additional 
damage  to  the  segment  of  the  armor  layer  damaged  by  Waves  4-6.  This  suggests  that 
the  damaged  armor  layer  can  withstand  subsequent  wave  action  at  a  lower  intensity  if 
the  wave  action  is  confined  to  the  damaged  area. 

As  one  would  expect,  the  damage  rates  differed  for  the  different  storm 
ordering.  But  the  ultimate  damage  at  the  end  of  the  storm  sequence  was  surprisingly 
consistent.  Both  the  mean  damage  and  the  standard  deviation  for  the  four  repeats  were 
consistent.  At  the  end  of  Series  B'  and  C  after  approximately  18,000  waves,  S=  8  and 
05  =  2  for  both  series.  Although  the  testing  was  by  no  means  comprehensive,  it  suggests 
that  damage  and  variability  are  similar  after  being  exposed  to  similar  cumulative  wave 
action  with  different  sequences.  This  is  useful  in  determining  the  reliability  of  a 
structure  because  the  cumulative  damage  may  be  assumed  to  be  the  sum  of  damages 
caused  by  individual  storms. 

5.5  Characteristics  of  Profile  Erosion 

For  Series  A',  the  mean  minimum  cover  depth  was  computed  as  the  average 
of  the  minima  from  the  16  profiles,  while  the  mean  eroded  depth  was  computed  as  the 
mean  of  the  maxima  from  the  16  profiles.  Figure  5.13  shows  the  mean  eroded  depth 
versus  number  of  waves,  and  Figure  5.14  shows  the  mean  cover  depth  versus  number  of 
waves,  both  for  Series  A'.  These  figures  clearly  show  the  periodic  healing  process  as  a 
jagged  sawtooth  shape  to  the  curves.  They  also  show  that  the  armor  layer  thickness 
decayed  to  an  equilibrium  level,  with  an  asymptote  at  C  =  0.  Figure  5.15  shows  that  the 
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eroded  length  jumps  to  a  level  of  L~1  initially  and  then  increases  with  the  number  of 
waves  in  a  manner  similar  to  damage,  as  one  would  expect.  The  initially  large  value  of 
L  computed  from  a  small  value  of  E  illustrates  that  there  was  quite  a  lot  of  stone 
settlement  but  not  displacement  as  damage  was  initiated. 
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Figure  5. 13.  Number  of  waves  versus  mean  maximum  eroded  depth  ±  one  standard 
deviation  for  Series  A ' 
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Figure  5. 14.  Number  of  waves  versus  mean  minimum  cover  depth  ±  one  standard 
deviation  for  Series  A' 
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Figure  5.15.  Number  of  waves  versus  mean  maximum  eroded  length  for  Series  A ' 
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Figure  5.16  shows  mean  eroded  depth  and  mean  cover  depth  versus  mean 
damage  for  Series  A'.  As  one  would  expect,  cover  depth  is  inversely  related  to  damage, 
while  eroded  depth  is  directly  related  to  damage.  Eroded  depth  and  cover  depth  are 
clearly  inversely  related.  Similarly,  Figure  5.17  shows  that  mean  eroded  length  and 
mean  damage  are  directly  related. 
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Figure  5.17.  Mean  damage  versus  mean  maximum  eroded  length  for  Series  A' 

Figures  5.18  through  5.23  show  that  the  profile  development  characteristics 
for  Series  B'mdC  are  similar  to  those  of  Series  A '  discussed  above.  Figures  5.24 
through  5.27  show  means  of  the  three  profile  parameters  versus  mean  damage  for  Series 
B’  and  C.  It  is  clear  that  profile  development  is  related  to  mean  damage.  This  will  be 
investigated  in  the  following  chapter. 
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Figure  5.18.  Number  of  waves  versus  mean  maximum  eroded  depth  ±  one  standard 
deviation  for  Series  B' 
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Figure  5.19.  Number  of  waves  versus  mean  maximum  eroded  depth  ±  one  standard 
deviation  for  Series  C 
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Figure  5.20.  Number  of  waves  versus  mean  minimum  cover  depth  ±  one  standard 
deviation  for  Series  B' 


Figure  5.21.  Number  of  waves  versus  mean  minimum  cover  depth  ±  one  standard 
deviation  for  Series  C 
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Figure  5.22.  Number  of  waves  versus  mean  maximum  eroded  length  for  Series  B' 


Figure  5.23.  Number  of  waves  versus  mean  maximum  eroded  length  for  Series  C 
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Figure  5.24.  Mean  damage  versus  mean  maximum  eroded  depth  and  mean  minimum 
cover  depth  for  Series  B' 


Figure  5.25.  Mean  damage  versus  mean  maximum  eroded  depth  and  mean  minimum 
cover  depth  for  Series  C 
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Figure  5.26.  Mean  damage  versus  mean  maximum  eroded  length  for  Series  B' 


Figure  5.27.  Mean  damage  versus  mean  maximum  eroded  length  for  Series  C 
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Chapter  6 

DAMAGE  AND  ERODED  PROFILE  PREDICTION 
6.1  Damage  Variability 

Figures  5.10  -  5.12  indicated  that  the  standard  deviation  of  damage  S 
varying  along  the  structure  increases  with  mean  damage  S.  Figure  6.1  shows  the 
relationship  between  O5  and  Sfor  all  data  in  the  three  series  plotted  in  Figures  5.10  - 
5.12.  The  scattered  data  points  for  the  three  series  may  be  represented  by  the  relation 

,T,=0is’“ 

which  is  the  solid  curve  shown  in  Figure  6.1.  This  relation  implies  that  damage 
variability  along  the  structure  increases  with  the  mean  damage.  On  the  other  hand, 
relative  variability  defined  as  the  coefficient  of  variation  V^  =  Og/ 1^  decreases  with 
increasing  S.  This  implies  that  higher  damage  levels  may  be  estimated  with  smaller 
relative  errors.  The  correlation  coefficients  for  the  three  series  for  Equation  6.1  were  all 
r  =  0.99.  An  example  of  the  use  of  Equation  6. 1  is  illustrated  for  5  =  13  at  failure  of 
Series  A'.  Substituting  S=  13  into  Equation  6.1  yields  =  2.65.  The  approximate 
range  of  -2.7  <  S*  <  3  in  Equation  5.10  corresponds  to  6  <  5  <  21,  indicating  damage 
variability  of  a  factor  of  3.5  along  the  structure. 
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•  Series  A'  ■  Series  B'  ▼  Series  C - Eq.  6.1 

> 

Figure  6.1.  Prediction  of  damage  variability,  characterized  by  the  standard  deviation,  as 
a  function  of  mean  damage  for  Series  A',  B',  and  C 

6.2  Eroded  Profile  Prediction 

It  is  useful  to  express  the  normalized  eroded  depth  Eand  the  variability  of 
eroded  depth,  a^,  as  a  function  of  mean  damage  also.  The  measured  values  of  and  E± 
Og  were  shown  as  a  function  of  cumulative  number  of  waves  for  Series  A '  in  Figure 
5.13.  The  trends  for  Series  B'  and  C  were  shown  to  be  similar  in  Figures  5.18  and  5.19. 
Figures  6.2  and  6.3  show  .^and  o^  as  a  function  of  Sfor  all  data  points  of  the  three 
series.  These  figures  indicate  the  following  empirical  relationships 


E  =0.46  VT 


(6.2) 


C7^  =  0.26-0.00007(5-7.8)"^ 


(6.3) 
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Figure  6.2.  Prediction  of  mean  maximum  eroded  depth  as  a  function  of  mean  damage 
for  Series  and  C 
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Figure  6.3.  Prediction  of  standard  deviation  of  maximum  eroded  depth  as  a  function  of 
mean  damage  for  Series  A',  B’,  and  C 


The  correlation  coefficients  for  Equation  6.2  were  r  =  0.99  for  all  three  series.  Equation 
6.2  indicates  EIS^'^  is  approximately  constant.  This  implies  that  the  average  shape  of 
the  eroded  area  remains  approximately  constant  during  damage  progression.  Figure  6.3 
and  Equation  6.3  suggest  that  the  variability  of  eroded  depth  increases  rapidly  with 
damage  and  becomes  approximately  constant  at  4.  The  variability  in  eroded  depth 
decreases  somewhat  as  the  damage  approaches  the  failure  level  of  the  armor  layer. 

The  eroded  length  L  along  the  slope  depicted  in  Figures  5.15, 5.17,  5.22, 
5.23,  5.26,  and  5.27  is  analyzed  in  the  following.  The  mean  eroded  length  4  is  simply 
defined  as  ^  =  2^,  /4>  corresponding  to  a  triangular  shape.  The  normalized  eroded 
length  L-Ig  is  then  given  by  Z  =  IS  IE,  as  shown  in  Equation  5.5.  From  the 
measured  values  of  ^and  E,  the  corresponding  values  of  Z  are  calculated  to  examine 
the  variation  of  Z  during  damage  progression.  Figure  6.4  shows  all  calculated  values  of 
Z  for  the  three  series  against  the  corresponding  values  of  S.  Figure  6.4  also  shows  the 
following  relationship  derived  from  Equations  6.2  and  5.5 

I  =4.4,,'X  <«■'» 

The  correlation  coefficients  for  Equation  6.4  were  r  =  0.98  for  Series  A'  and  r  =  0.99  for 
Series  B’  and  C.  As  an  example  of  the  use  of  Equation  6.4,  for  5  =  13  at  failure  for 
Series  A’,  Equations  6.2  and  6.4  yield  E~\.l  and  Z  =  16,  implying  the  average 
damaged  profile  at  failure  extends  16  stone  diameters  along  the  slope.  For  the  model 
structure  with  D„5o  =  3.64  cm  shown  in  Figures  4.2  and  4.3,  this  damage  extent 
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corresponds  to  most  of  the  seaward  slope  and  crest  of  the  structure  as  shown  in  Figures 
5.8  and  5.9. 


•  Series  A'  «  Series  B'  ▼  Series  C - Eq.  6.4 

Figure  6.4.  Prediction  of  mean  eroded  length  as  a  function  of  mean  damage  for  Series 
A',  5',  and  C' 

Finally,  the  normalized  minimum  cover  depth  C  is  analyzed  to  examine  the 
decrease  of  protection  with  damage  progression.  The  measured  values  of  C  and  C± 
were  plotted  as  a  function  of  cumulative  number  of  waves  in  Figures  5.14, 5.20,  and 
5.21  for  Series  A',  B’,  and  C,  respectively.  These  relations  are  similar  to  the 
corresponding  results  for  damage  and  eroded  depth  except  that  C  decreases  with  N^.  and 
Oc  is  positive  at  N„  =  0.  The  initial  values  of  C  at  A/,,,  =  0  were  Q  =1.65, 1.52,  and  1.38 
for  Series  A',  B’,  and  C',  respectively,  and  the  initial  value  of  Oc  at  =  0  was  Oco  = 

0. 14, 0. 17,  and  0.22  for  Series  A’,B',  and  C,  respectively.  To  account  for  the 
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differences  of  Q  and  a^o,  (Q  -  Q  and  (Oc  -  0^)  are  plotted  against  i^for  all  the  data 
points  of  the  three  series  in  Figures  6.5  and  6.6,  respectively.  The  empirical 
relationships  shown  in  these  figures  are  given  by 


C  =  Co -0.15 

<^c  =^co+0.O98 -0.002(5  -7)^ 


(6.5) 

(6.6) 
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•  Series  A'  ■  Series  B'  ▼  Series  C - Eq.  6.5 

Figure  6.5.  Prediction  of  mean  minimum  cover  depth  as  a  function  of  mean  damage  for 
Series  A',  5',  and  C 
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•  Series  A'  ■  Series  B'  ▼  Series  C - Eq.  6.6 


Figure  6.6.  Prediction  of  standard  deviation  of  minimum  cover  depth  as  a  function  of 
mean  damage  for  Series  A',  B',  and  C 

For  Equation  6.5,  the  correlation  coefficients  were  r  =  0.99  for  all  three 
series.  Figure  6.5  and  Equation  6.5  indicate  that  the  mean  cover  depth  decreases 
approximately  linearly  with  mean  damage  increase.  The  scattered  values  of  (Oc  - 
plotted  in  Figure  6.6  are  of  the  same  order  of  magnitude  as  the  initial  value  Oco  listed 
above.  Figure  6.6  and  Equation  6.6  indicate  that  the  standard  deviation  representing 
the  variability  of  cover  depth  along  the  structure  increases  somewhat  with  damage 
progression  and  then  decreases  as  the  cover  depth  approaches  zero.  It  is  noted  that 
Equation  6.6  overestimates  o^  somewhat  for  Series  C  as  can  be  seen  in  Figure  6.6. 

A  failure  criterion  based  on  the  normalized  minimum  cover  depth  C  is 
explored  in  the  following.  Using  the  lower  limit  of  C*  =  (C  -  Q  lOc  in  Equation  5.10, 
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the  lower  limit  of  C  may  be  estimated  as  (C  The  failure  may  be  assumed  to 

initiate  when  this  lower  limit  becomes  zero.  This  criterion  with  Equations  6.5  and  6.6 
yields  the  following  equation  for  the  mean  damage  at  failure  initiation 

0.15  +  2.7[0.098  -  0.002(5-7)2]  =  ^  -  2.1  a (6.7 

This  criterion  accounts  for  the  initial  minimum  cover  depth  and  its  variability  along  the 
structure  where  (Q  -2.70co)  =  1.27, 1.06,  and  0.79  for  Series  A',  B',  and  C.  Solving 
Equation  6.7  for  mean  damage  for  Series  A'  yields  5=  10.9  at  failure  initiation, 
compared  with  5=  13  measured.  The  failure  initiation  criterion  of  Equation  6.7  based 
on  the  lower  limit  of  C*  of  all  data  points  plotted  in  Figure  5.4  yields  the  lower  limit  of 
the  mean  damage  at  failure.  The  upper  limit  of  mean  damage  at  failure  may  simply  be 
estimated  by  C  =  0  in  Equation  6.5.  This  upper  limit  is  5=  16.5  for  Series  A'. 

6.3  Temporal  Damage  Development 

If  the  mean  damage  can  be  predicted,  then  Equations  6.1-  6.6  yield  Oj,  E, 
Of,  L,  C,  and  o^.  These  relationships  have  been  obtained  by  analyzing  the  measured 
profiles  of  the  armor  layer  and  underlayer  without  regard  to  the  incident  waves  and 
water  level.  This  statistical  analysis  indicates  that  the  damaged  profile  statistics  can  be 
represented  by  the  mean  damage  5.  The  next  step  is  to  predict  the  temporal  variation  of 
5  on  the  rubble-mound  breakwater  exposed  to  depth-limited  breaking  waves  in 
sequences  of  storms  with  varying  wave  conditions  and  water  levels. 
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The  empirical  formulas  of  Van  der  Meer  (1988)  indicate  that  breakwater 
damage  due  to  depth-limited  waves  is  a  function  of  N^,  T„,  and  time  t.  In  addition,  it  is 
likely  that  damage  is  a  function  of  structure  slope,  structure  permeability,  armor 
gradation,  armor  porosity,  armor  stone  shape,  and  method  of  armor  placement.  None  of 
these  latter  parameters  except  armor  gradation  were  varied  in  this  study.  Gradation  is 
discussed  in  Chapter  7.  The  remainder  of  the  parameters  will  be  investigated  in  future 
studies. 

In  the  following,  the  depth-limited  breaking- wave  stability  data  of  Van  der 
Meer  (1988)  will  be  used  to  develop  relations  for  damage  progression.  His  depth- 
limited  wave  and  measured  damage  data  are  shown  in  Table  6.1.  The  particular  data  set 
shown  in  Table  6.1  corresponded  to  a  structure  slope  of  1V:2H,  a  beach  slope  of 
1V:30H,  and  the  stone  armor  characterized  by  Dgs  /  £>15=  1.25,  I>„5o=  3.60  cm,  and 
A=1.615.  The  table  shows  that,  for  toe  depth  h,  =  0.4  m,  the  waves  were  only 
marginally  depth  limited  because  H^lh,  was  in  the  range  0.26-0.3 1 .  Therefore,  Van  der 
Meer’s  data  provide  only  eight  tests  where  the  waves  were  clearly  depth  limited  Qi,  = 

0.2  m).  The  temporal  development  of  damage  cannot  be  deduced  from  these  data 
because  the  damage  was  only  given  after  1,000  waves  and  3,000  waves. 

From  Van  der  Meer’s  data,  it  can  be  seen  that  the  surf  similarity  parameter 
with  respect  to  the  beach  slope  for  the  eight  breaking  wave  tests  9-16  ranged  from  0.20 
to  0.25.  This  suggests  that  all  were  spilling  breakers,  which  is  not  the  worst  case  for 
armor  stability.  As  explained  in  the  previous  chapters,  plunging  breakers  expose  the 
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armor  units  to  the  highest  vertical  forces  at  the  steep  front  of  a  plunging  wave,  which  is 
the  worst  case  for  armor  stability. 


Table 

Depth 

6.1 

i-Limited  Wave  Damage  Data  from  Van  der  Meer 

f1988) 

Test 

Toe 

Depth 

m 

Wave 

Height 

m 

Mean 

Per¬ 

iod 

T„, 

sec 

HJh, 

Ns 

^tor 

Beach 

Slope 

^tor 

Struc¬ 

ture 

Slope 

Sfor 

1,000 

1 

0.4 

0.1253 

1.68 

0.31 

2.16 

0.198 

2.97 

8.51 

12.69 

2 

0.4 

0.0900 

1.69 

0.23 

1.55 

0.235 

3.52 

0.86 

1.56 

3 

0.4 

0.1088 

1.68 

0.27 

1.87 

0.212 

3.18 

2.85 

4.05 

■1 

0.4 

0.1359 

1.67 

0.34 

2.34 

0.189 

2.83 

10.15 

18.15 

5 

0.4 

0.1286 

2.19 

0.32 

2.21 

0.255 

3.82 

6.04 

9.85 

6 

0.4 

0.1052 

2.19 

0.26 

1.81 

0.281 

4.22 

1.16 

2.15 

■ 

0.4 

0.1555 

2.18 

0.39 

2.67 

0.231 

3.46 

8.78 

20.04 

8 

0.4 

0.1215 

2.20 

0.30 

2.09 

0.263 

3.94 

4.78 

7.19 

9 

m 

0.1345 

2.18 

0.67 

2.31 

0.248 

3.72 

4.50 

6.29 

10 

0.2 

0.1275 

2.17 

0.64 

2.19 

0.253 

3.80 

3.05 

3.00 

11 

0.2 

0.1415 

2.17 

0.71 

2.43 

0.241 

3.61 

5.13 

6.43 

12 

0.2 

0.1445 

2.18 

0.72 

2.49 

0.239 

3.58 

6.52 

12.45 

13 

0.2 

0.1270 

1.70 

0.64 

2.18 

0.199 

2.98 

2.27 

3.63 

14 

0.2 

0.1320 

1.73 

0.66 

2.27 

0.199 

2.98 

4.97 

7.97 

15 

0.2 

0.1350 

1.74 

0.68 

2.32 

0.197 

2.96 

6.65 

11.03 

16 

0.2 

0.1350 

1.77 

0.68 

2.32 

0.201 

3.01 

8.65 

11.38 

Consider  Equations  3.7  and  3.8,  which  are  Van  der  Meer’s  equations  for 
stability  that  include  damage  as  a  parameter.  The  stability  number  was  a  minimum  at 
the  transition  between  plunging  and  surging  waves  occurring  at 
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=  [6.2P®-^’v/tan0]2 


(6.8) 


where  g  =  (^  +  0-5)''  and  P  =  permeability  coefficient  varying  in  the  range  0. 1  -  0.6. 
Keep  in  mind  that  0  in  Equation  6.8  is  the  structure  slope  measured  from  horizontal. 
Later  in  this  chapter,  this  requirement  will  be  abandoned.  Using  Equations  3.7  and  6.8, 
the  minimum  stability  number  can  be  expressed  as 


N. 


=  p(2(O.18/’-O.O65)^j,Qj0)O.25e 


(6.9) 


where 

N,  =  H,/AD„5o  =  stability  number  based  on  the  significant  wave  height 
C;v  =  empirical  coefficient  introduced  here 

Ny,  =  number  of  waves  associated  with  the  damage  S  starting  from  zero 
damage 


For  depth-limited  breaking  waves  on  a  sloping  beach  in  front  of  the 
structure.  Van  der  Meer  (1988)  proposed  the  use  of  H2%  and  rewrote  his  equation  for 
relatively  deep  water  using  H2%  =  based  on  the  Rayleigh  distribution  of  wave 
heights  where  based  on  H,  was  not  changed.  Note  that  the  values  of  H2%  were  not 
tabulated  in  his  report,  so  it  is  difficult  to  evaluate  this  shallow-water  extension  of  his 
stability  equations.  In  the  present  experiment,  incident  irregular  wave  breakers  were 
spilling  and  plunging  on  the  1:20  beach  slope  and  collapsing  and  surging  on  the  1:2 
structure  slope.  However,  plunging  waves  on  the  beach  were  also  observed  to  hit  the 
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structure  directly.  Consequently,  although  it  is  likely  that  the  beach  slope  was  most 
influential  on  the  breaking  process,  an  appropriate  slope  for  is  uncertain  and  Equation 
6.9  is  simply  adopted,  although  the  relation  does  not  contain  The  stability  number 

N,  in  Equation  6.9  based  on  is  used  here  because  it  will  be  easier  to  predict  than 
i/2%  using  time-averaged  surf  zone  models  such  as  Battjes  and  Stive  (1985)  and  Dally 
(1992)  that  will  need  to  be  modified  to  account  for  reflected  waves  (Baquerizo  et  al. 
1997).  The  empirical  coefficient  is  calibrated  using  the  16  tests  in  Table  6.1.  For 
eight  tests  with  toe  depth  h,  =  0.4  m  and  /h,  =  0.23-0.39,  Q  =  0.97-1.21  for  N„  = 

1 ,000  and  =  0.99- 1 . 19  for  =  3,000.  For  eight  tests  with  h,  =  0.2  m  and  /h,  = 

O. 64-0.72,  Cf/  =  1.03-1.27  for  =  1,000  and  =  1.09-1.34  for  N,.  =  3,000.  The 

present  experiment  with  H,  th,  =  0.64-1.1 1  using  Table  4.4  is  more  similar  to  the  eight 
tests  in  Table  6. 1  with  h,  =  0.2  m.  So  Cf,=  \.2  is  tentatively  assumed  in  Equation  6.9. 

Equation  6.9  is  rewritten  as  two  damage  formulas  in  the  following  to 
facilitate  the  representation  of  incident  random  waves  in  time  series  and  spectra.  The 
first  damage  model,  based  on  wave  time  series  statistics,  has  the  form 


where  t=T^  is  the  test  duration  for  constant  wave  conditions.  The  empirical 
coefficient  a,  is  related  to  cot  0,  P,  and  C„,  and  b  in  Equation  6. 10  is  introduced  for  long 
duration  tests,  where  b  =  0.5  in  Equation  6.9.  The  second  model  is  based  on  the  wave 
frequency  domain  statistics  and  has  the  form 
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(6.11) 


In  Equation  6.1 1,  /AD„5o  is  the  stability  number  based  on  the  zeroth  moment 

wave  height,  the  test  duration  t  for  constant  wave  conditions  is  normalized  by  the  peak 
period,  and  and  b  are  again  empirical  coefficients  that  will  be  obtained  from  data.  It 
is  noted  that  Equations  6.10  and  6.1 1  assume  5  =  0  at  f  =  0  and  constant  wave 
conditions. 


The  expression  for  in  Equation  6.10  can  be  found  from  Equation  6.9 
where  P  =  0.4  for  conventional  rubble-mound  breakwaters.  For  Q  =  1.2,  P  =  0.4,  and 
cot  0  =  2,  Uj  =  0.003  for  Equation  6.10  ifb  =  0.5  is  applicable  to  the  present  long 
duration  tests.  Van  der  Meer  (1988)  analyzed  the  five  long  duration  tests  with  N^.  up  to 
15,000  conducted  by  Thompson  and  Shuttler  (1976)  and  obtained  the  term  S IN^-^  in 
Equation  6.9  for  =  1,000-8,500.  His  data  analysis  indicated  b  <  0.5  for  >  8,500. 
In  the  following,  the  values  of  a,  and  b  are  calibrated  for  this  experiment  where  Ny,= 
60,000  for  Series  A'  and  iV,,  =  18,000  for  Series  B'  and  C'. 

The  empirical  relationship  in  Equation  6.10  is  limited  to  damage  due  to 
incident  irregular  waves  with  constant  and  T„  starting  from  1?=  0  at  r  =  0.  Asa 
result,  these  formulas  and  other  existing  formulas  are  intended  for  prediction  of  damage 
during  the  peak  of  a  design  storm.  In  order  to  develop  an  empirical  procedure  for  more 
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realistic  conditions  of  and  varying  with  time,  the  rate  of  damage  increase  dSIdt  is 
obtained  from  Equation  6.10  as 

^  (6.12) 
at 

which  is  assumed  to  be  valid  at  arbitrary  time  t.  To  apply  Equation  6.12  for  and  T„ 
varying  with  time  t,  the  relation  must  be  integrated  numerically  with  an  arbitrary  initial 
value  of  S.  For  practical  applications,  the  values  of  and  may  be  assumed  to  be 
constant  for  a  short  duration.  Integration  of  Equation  6.12  for  the  duration  of  constant 
and  lasting  r  =  4  to  r  =  yields  the  mean  damage  at  arbitrary  time  t 

5(0  =  5(r„)  -  tf)  for  (6.13) 

where  S(t„)  =  known  damage  at  t  - 1„.  Use  of  Equation  6.13  for  each  interval  of 
constant  and  T„  in  sequence  allows  one  to  compute  Sit„)  for  incident  wave 
conditions  represented  by  H,  and  varying  with  time.  Alternatively,  representing  the 
irregular  waves  at  the  toe  of  the  slope  using  spectral  parameters  as  shown  in  Equation 
6.11  yields  the  equation 

5(0  =  5(t„)  -  a^NlT;\t »  - 1^)  for  <  t  <  (6.14) 

Equations  6.13  and  6.14  are  valid  for  an  undamaged  structure  (5  =  0  at  r  =  0)  or 
damaged  structure,  S{tf}  =  5^  at  r  =  t^.  Using  the  coefficients  predicted  by  Van  der 
Meer’s  data,  =  0.003  and  b  =  0.5,  Equation  6.13  follows  the  general  trend  of  the  data 
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but  underpredicts  damage  up  to  A^„~5000.  For  a  greater  number  of  waves,  this  equation 
overpredicts  damage. 

The  empirical  coefficients  in  Equations  6.13  and  6.14  are  recalculated  in  the 
following  using  the  present  experimental  data,  which  were  obtained  from  tests  much 
longer  in  duration  than  those  listed  in  Table  6.1.  The  wave  heights  and  periods  for  the 
present  experiment  were  listed  in  Tables  4.3  and  4.4  and  the  experimental  data  are 
summarized  in  Appendix  B.  In  order  to  get  a  more  gradual  increase  in  S(t)  it  was 
necessary  to  raise  a  and  lower  b.  The  final  coefficients  were  a,  =  0.025  andb  =  0.25  for 
Equation  6. 13  and  =  0.022  andb  =  0.25  for  Equation  6. 14  yielding  the  general 
predictive  relations  for  mean  damage  progression  due  to  breaking  waves  for  time  and 
frequency  wave  statistics,  respectively,  as 

Sit)  =  Sit„)  +  0.025(/V,)-'C for  (6.15) 

5(0  =  5(t„)  +  0.022(iV„/r/-""  for  t„<t^  (6.16) 

Figure  6.7  shows  Equation  6.15  (dashed  line)  and  Equation  6.16  (solid  line)  plotted 
against  the  damage  data  of  Figure  5. 10  for  Series  A '.  The  mean  damage  is  well 
predicted  by  Equations  6.15  and  6.16  with  a  correlation  coefficient  of  r  =  0.99  for  both 
equations  for  this  series.  Similarly,  Figures  6.8  and  6.9  show  that  the  equations  predict 
damage  well  for  Series  B'  and  C',  with  correlation  coefficients  of  r  =  0.99  for  both. 
Equations  6.15  and  6. 16  are  therefore  generalized  for  Series  A',  B',  and  C.  The 
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agreement  for  Equations  6.15  and  6.16  is  good  for  the  three  series  except  that  the 
formulas  do  not  predict  damage  stabilization  adequately.  This  deficiency  becomes 
much  worse  for  b  =  0.5.  The  calibrated  value  of  ^  =  0.25  indicates  that  damage  varies 
less  with  time  t  in  the  long  duration  experiments  using  uniform-sized  stone  armor.  The 
calibrated  values  of  =  0.025  and  Cp  =  0.022  for  b  =  0.25  are  much  larger  than  = 
0.003  for  b  =  0.5  estimated  using  data  from  Table  6. 1 .  This  difference  is  mainly  caused 
by  the  different  value  of  b.  As  an  example,  for  t  /T^  =  10,000,  (t  /Tj’  =  10  or  100  for 
b  =  0.25  or  0.5.  Consequently,  the  values  of  and  Op  are  coupled  with  the  calibrated 
value  of  b. 


Figure  6.7.  Damage  prediction  relations  compared  to  data  for  Series  A' 
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2.2 


2.2 


2.5 


Figure  6.8.  Damage  prediction  relations  compared  to  data  for  Series  B' 


Figure  6.9.  Damage  prediction  relations  compared  to  data  for  Series  C 
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Chapter  7 


EFFECT  OF  WAVE  PERIOD  AND 
ARMOR  GRADATION  ON  DAMAGE  PROGRESSION 


7.1  Experimental  Setup  and  Test  Conditions 

It  is  expected  that  the  empirical  coefficients  in  Equations  6.15  and  6.16  are, 
to  some  extent,  a  function  of  wave  period.  In  addition,  only  uniform  sized  armor  was 
utilized  in  Series  A',  B\  and  C’.  Therefore,  four  test  series  (£)',  E',  F',  and  G")  were 
added  to  address  the  effect  of  wave  period  and  stone  gradation  on  damage  progression. 
These  additional  test  series  are  summarized  in  Table  7.1  along  with  Series  A',  B',  and  C 
for  comparison.  The  spectral  and  time  series  incident  wave  parameters  for  all  series  are 
listed  in  Table  7.2.  Damage  measurements  are  summarized  for  Series  D\ E\  F\  and  G' 
in  Tables  B4  through  B7  in  Appendix  B. 


1  Table  7.1 

I  Summary  of  Test  Series  for  Damage  Ex 

periments 

Test 

Series 

Test 

Type 

Armor 

Type 

Water  Level 
Order 

Test  Duration 
(hr) 

A' 

Deterioration  to  Failure 

Uniform 

Low  -  High 

28.5 

B' 

Storm  Ordering 

Uniform 

Low  -  High 

8.5 

C' 

Storm  Ordering 

Uniform 

High  -  Low 

9.0 

D' 

Wave  Period 

Uniform 

Low-High 

8.5 

E' 

Wave  Period 

Uniform 

Low-High 

8.5 

mm 

Gradation 

Riprap 

Low-High 

8.5 

G' 

Gradation 

Riorao 

8.5 
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Table  7.2 

Summary  of  incident  Wave  Characteristics 


erie 

Wave 

Dura¬ 

tion 

hr 

Ms 

Ml/to 

Mgx 

cm 

cm 

cm 

1 

1.5 

11.9 

2.48 

9.78 

9.38 

■  11.50 

2 

1.5 

11.9 

2.48 

12.40 

0.47 

1.69 

11.60 

13.80 

15.37 

3 

n 

11.9 

2.48 

14.20 

0.48 

1.74 

13.20 

15.70 

17.30 

4 

1.0 

15.8 

2.59 

10.50 

0.52 

1.73 

10.10 

12.72 

14.30 

5 

15.8 

2.59 

13.60 

0.51 

1.67 

13.00 

15.97 

17.80 

6 

11.0 

15.8 

2.59 

15.80 

0.51 

1.66 

■El^ 

18.00 

19.30 

11.9 


15.8 


15.8 


1.0  15.8 


15.8 


15.8 


11.9 


15.80 


10.50 


13.60 


15.80 


12.40 


14.20 


6.13 


9.88 


13.11 


9.62 


12.83 


9.38 


11.60 


13.20 


13. 


14.90 


10.10 


13.00 


14.90 


11.60 


13.20 


6.05 


9.88 


13.18 


9.74 


13.21 


5.05 


7.26 


10.19 


11.50 


13.80 


15.70 


15.97 


18. 


12.72 


15.97 


18. 


13.80 


15.70 


12.80 


15.37 


17.30 


17.80 


19.30 


14.30 


17.80 


19.30 


15.37 


17.30 


8.36 


13.59 


17.11 


14. 


17.87 


7.75 


11.11 


14.90 


8.78 


13.34 


9.61 


15.63 


19.39 


8.58 


12.45 


10.62 


16.22 


18.99 


16.67 


20.03 


11.9 

1.53 

5.05 

0.38 

1.29 

11.9 

1.53 

7.13 

0.35 

1.29 

11.9 

1.53 

9.93 

0.31 

1.23 

15.8 

1.48 

6.60 

0.34 

1.30 

15.8 

1.48 

9.41 

0.32 

1.26 

12.48 


16.14 


12.48 


16.80 


6.72 


9.70 


13.38 


8.15 


12.03 


8.68 


14.39 


18.09 


7.78 


10.95 


9.53 


14.93 


17.95 


15.06 


18.64 


9.53 


6.96 


11.51 


14.95 


6.18 


8.54 


7.53 


11.99 


15.21 


7.21 


11.68 


15.33 


6.43 


8.82 


11.9 

1.97 

7.62 

0.42 

1.50 

11.9 

1.97 

12.07 

0.37 

1.36 

11.9 

1.97 

15.42 

0.35 

1.30 

15.8 

2.02 

11.92 

0.37 

1.44 

15.8 

2.02 

15.34 

0.35 

1.34 

15.36 


Series  D'  and  E'  were  similar  to  Series  B'  except  with  different  target  peak 
periods.  The  different  spectra  also  produced  different  characteristic  wave  heights  at  the 
toe  of  the  structure.  Series  F’  and  G’  were  similar  to  Series  B'  and  D'  except  that  riprap 
was  used  as  armoring.  Series  F'  had  the  same  wave  peak  periods  as  Series  B’  while 
Series  G'  had  the  same  wave  peak  periods  as  Series  D'.  Therefore,  Series  B',  D',  E',  F', 
and  G'  provide  comparative  damage  measurements  with  systematic  variation  of  armor 
gradation  and  wave  period. 

As  shown  in  Table  7.1,  the  order  of  water  levels,  individual  storm  durations, 
and  series  durations  for  Series  D',  E',  F\  and  G'  were  identical  to  Series  B'.  Four  storms 
were  run  within  each  test  series  following  30  min  of  shakedown  waves,  with  the  low- 
water-level  storms  run  first  followed  by  the  high-water-level  ones.  The  duration  of  each 
of  4  storms  was  2  hours  yielding  8.5  hours  of  waves  for  each  series.  Wave  heights  were 
incrementally  raised  at  each  water  level  to  simulate  increasing  storm  severity.  The 
waves  were  run  in  15  min  bursts,  with  the  water  completely  settling  between  bursts.  As 
was  done  for  Series  A',  B',  and  C',  the  2  structures  were  profiled  every  30  min  using 
eight  profile  rods  across  the  center  portion  of  each  structure.  Each  storm  series  was 
repeated  to  yield  32  profiles  alongshore  for  the  undamaged  underlayer  and  armor  layer 
and  after  every  30  min  of  waves.  As  discussed  in  Chapter  5,  the  mean  and  standard 
deviation  of  E,  C,  L,  and  5  for  these  32  profiles  were  computed. 

The  armor  stone  for  Series  F'  and  G'  was  widely  graded  riprap  with  a 
median  mass  Mjo  =  256  g,  nonrinal  diameter  =  4.58  cm,  stone  density 
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=  2.66  g/cm^  and  =  1.53.  The  mass  distribution,  plotted  as  percent  finer  by 

mass,  is  shown  in  Figure  7. 1 .  The  median  mass  was  computed  from  a  200  stone  sample 
as  the  value  where  half  the  total  mass  in  the  sample  was  greater.  For  comparison,  the 
median  mass  by  count  for  a  200  stone  sample  was  =  128  g  (i.e.  half  the  stones 
weighed  less).  This  latter  value  was  not  used  in  this  thesis,  although  it  is  sometimes 
used  in  the  field.  The  riprap  followed  the  widest  recommendation  of  the  SPM  (1984)  of 
approximately  0.125M5o  <  M  <  For  all  series,  the  underlayer  or  filter  layer  had  a 
gradation  of  =  1.32  and  was  sized  such  that  /  (Afso)^;,,^,  =  25  and 

(Ps^armor  !  P^filter  =  2.9. 


Figure  7.1.  Stone  mass  distributions  for  riprap  and  underlayer 
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7.2  Probability  Density  Functions  of  Damage,  Eroded  Depth,  and  Cover  Depth 

The  probability  density  functions  of  S*,  E*,  and  C*  are  shown  in  Figures 
7.2, 7.3,  and  7.4,  respectively,  with  Series  D’,  E',  F'  and  G'  combined  onto  each  figure. 
These  figures  can  be  compared  to  Figures  5.2, 5.3,  and  5.4  for  Series  A',  B',  and  C', 
respectively.  The  figures  indicate  that  the  variability  of  the  eroded  profile  and  eroded 
area  for  the  new  series  are  similar  to  those  of  Series  A',  B',  and  C.  The  limits  of 
variability  remain  as  prescribed  by  Equation  5.10  or  approximately  in  the  range  from 
-3.0  to  3.0.  It  is  noted  that  the  damage  and  eroded  depth  distributions  are  negatively 
skewed  while  the  cover  depth  is  positively  skewed.  This  indicates  that  more  cross 
sections  have  less  normalized  damage  and  larger  normalized  cover  layer  thickness. 
Relatedly,  the  negative  tails  of  the  damage  and  eroded  depth  distributions  are  less 
populated  for  Series  D\  E',  F'  and  G'  than  for  Series  A',  B',  and  C',  particularly  below 
normalized  values  of  -2.  As  will  be  described  in  the  following  sections,  wave  energy 
and  therefore  damage  levels  were  lower  for  the  new  series  than  for  Series  A',  B',  and  C'. 
The  damage  increased  rapidly  at  low  damage  levels,  5=0-1,  resulting  in  few  damage 
measurements  at  these  low  damage  levels. 
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Figure  7.2.  Probability  density  function  for  normalized  damage  S*  for  Series  D',  E\  F' 
andG' 


Figure  7.3.  Probability  density  function  for  normalized  eroded  depth  E*  for  Series  D\ 
F,F'andG' 


132 


c* 


•  Series  D'  “  Series  E'  ■»  Series  F'  »  Series  G' - Normal 


Figure  7.4.  Probability  density  function  for  normalized  cover  depth  C*  for  Series  D',  E\ 
F'andG' 


7.3  Damage  Variability 

Figure  7.5  shows  Og  as  a  function  of  mean  damage  5 for  Series  D\  E',  F\ 
and  G'  along  with  Equation  6.1.  Series  B’  is  included  for  comparison.  This  figure 
indicates  that  Equation  6. 1  underpredicts  damage  variability  for  the  new  series  in  the 
range  S=  1-  4.  But  above  and  below  this  range  the  results  are  inconclusive.  The  greater 
variability  in  damage  for  the  wider  gradation  (Series  F’  and  G”)  is  expected;  but  the 
reason  for  the  greater  damage  variability  for  shorter  wave  periods  is  not  clear.  It  appears 
there  are  insufficient  data  to  make  a  definitive  modification  to  Equation  6.1.  Even 
though  the  data  are  slightly  underpredicted  by  Equation  6.1,  Equations  5.10  and  6.1 
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provide  a  description  of  the  alongshore  variability  of  damage  that  is  reasonable  for 
experimental  results  described  herein. 


•  Series  B'  ■  Series  D’  ▼  Series  E’ 

*  Series  F'  *■  Series  G' - Eq.  6.1 


Figure  7.5.  Prediction  of  damage  variability,  characterized  by  the  standard  deviation,  as 
a  function  of  mean  damage  for  Series  B',  D',  E\  F',  and  G' 

7.4  Eroded  Profile  Prediction 

Figure  7.6  shows  the  mean  eroded  depth  J^for  the  four  new  series  as  a 
function  of  5^ along  with  Series  B’.  Equation  6.2,  for  prediction  of  E  based  on  Series  A', 
B\  and  C,  is  also  plotted  in  Figure  7.6.  It  is  clear  that  this  equation  describes  the  new 
data  well  without  modification.  The  standard  deviation  for  eroded  depth  as  a 
function  of  i'is  shown  in  Figure  7.7  for  the  four  new  series  and  B’.  Equation  6.3  for 
as  a  function  of  S,  also  plotted  in  Figure  7.7,  underestimates  the  variability  of  eroded 


134 


depth  in  the  new  data  sets.  This  reinforces  the  fact  that  damage  variability  was  greater 
in  Series  £)',  E',  F',  and  G'  than  in  Series  A',  B’,  and  C.  Also,  Series  D'  produced  more 
outlying  points  in  this  Figure  than  the  other  series.  This  is  consistent  with  previous 
figures  and  indicates  that  further  testing  is  required  to  better  quantify  the  coefficient  of 
variation  of  the  damage  and  profile  parameters,  particularly  for  low  damage  levels.  The 
mean  normalized  eroded  length  L  is  plotted  as  a  function  of  .y  for  data  from  the  four 
new  series  in  Figure  7.8.  Equation  6.4  for  Z  as  a  function  of  i^is  shown  to  provide  an 
excellent  fit  to  the  new  data. 


•  Series  B'  «  Series  D'  ▼  Series  E' 
+  Series  F'  *  Series  G' - Eq.  6.2 


Figure  7.6.  Prediction  of  mean  maximum  eroded  depth  as  a  function  of  mean  damage 
for  Series  B',  £>',  £',  F',  and  G’ 
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•  Series  B'  "  Series  D'  »  Series  E' 

+  Series  F'  *  Series  G’ - Eq.  6.3 

Figure  7.7.  Prediction  of  standard  deviation  of  maximum  eroded  depth  as  a  function  of 
mean  damage  for  Series  B\ D',  E',  F\  and  G' 


Mean  Damage,  S 

•  Series  B'  ■  Series  D'  ▼  Series  E’ 

+  Series  F'  a  Series  G' - Eq.  6.4 

Figure  7.8.  Prediction  of  mean  maximum  eroded  length  as  a  function  of  mean  damage 
for  Series  B',D',  E',  F',  and  G' 
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The  mean  cover  depth  as  a  function  of  5  was  shown  to  be  described  by 
Equation  6.5.  This  relation,  rearranged  to  be  eroded  cover  depth,  also  provides  a  good 
fit  to  the  new  data,  as  shown  in  Figure  7.9.  In  addition,  the  standard  deviation  of 
minimum  cover  depth  Oc  was  shown  to  be  described  by  Equation  6.6.  This  equation, 
along  with  the  new  data,  are  shown  in  Figure  7.10.  Equation  6.6  underpredicts  the 
variability  in  cover  depth,  as  expected.  Again,  Series  D'  produced  most  of  the  outlying 
points.  Data  from  the  other  series  are  reasonably  well  predicted  by  Equation  6.6. 

The  profile  and  variability  relationships  given  by  Equations  5.10  and  6.1 
through  6.6,  derived  from  Series  A',  B',  and  C  data,  fit  the  new  data  remarkably  well, 
considering  the  ranges  of  stone  grading  and  peak  wave  period  in  this  experiment. 
Equation  6.3  and  6.6  underpredict  the  variability  and  will  likely  require  modification  as 
more  data  are  acquired.  It  appears  that  the  coefficient  of  variation  is  more  scattered  for 
low  damage  levels.  The  relations  given  above  allow  prediction  of  profile  shape  and 
alongshore  variability  of  the  profile.  The  above  relations  were  developed  using  only 
profile  data  without  regard  to  incident  wave  conditions  or  water  levels. 
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Figure  7.9.  Prediction  of  mean  minimum  cover  depth  as  a  function  of  mean  damage  for 
Series  E',  F\  and  G' 


Figure  7.10.  Prediction  of  standard  deviation  of  minimum  cover  depth  as  a  function  of 
mean  damage  for  Series  B’,  D\  E',  F’,  and  G' 
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7.5  Temporal  Damage  Development 

In  Chapter  6,  the  generalized  Equations  6.15  and  6.16  for  predicting  damage 
progression  were  highly  correlated  with  damage  data  from  Series  A',  B',  and  C'.  Figures 
7.1 1  through  7.14  show  Equations  6.15  and  6.16  plotted  against  data  from  Series  D',  E\ 
F',  and  G'.  To  recap,  Series  D'  and  E'  were  nearly  identical  to  Series  B'  except  that  the 
peak  wave  period  was  changed  for  the  three  series.  Series  F'  and  G'  were  again  similar 
except  that  the  uniform  armor  was  replaced  with  riprap.  Two  different  peak  wave 
periods  were  tested  in  Series  F'  and  G'.  It  can  be  seen  that  the  damage  progression 
equations  predict  the  overall  trend  of  damage  progression  for  Series  D'  through  G'  using 
=  0.025,  Up  =  0.022,  and  b  =  0.25,  although  there  are  noted  discrepancies.  For 
example,  it  can  be  seen  that  damage  initiation  is  underpredicted  in  all  series.  Equations 
6.15  and  6.16  significantly  underpredict  damage  initiation,  if  only  1  or  2  stones  are 
displaced  at  the  beginning  of  each  test  series.  This  underprediction  appears  to  be 
produced  by  the  variability  in  damage  initiation.  Figures  7.15  through  7.18  show 
Equation  6. 15  starting  from  two  different  points:  from  =  0  and  from  the  first 
measured  damage  point  past  N,,.  =  1000.  As  can  be  seen,  the  prediction  is  much  better 
when  the  damage  is  predicted  after  initial  profile  adjustment.  Figures  7.1 1  through  7.18 
indicate  that  the  empirical  coefficients  in  Equations  6.15  and  6.16  may  vary  somewhat 
with  wave  period  and  stone  gradation  for  the  range  of  experimental  conditions  described 
herein.  The  initial  profile  adjustment  may  need  to  be  accounted  for  in  test  series  with 
relatively  small  cumulative  damage. 
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Figure  7.11.  Damage  prediction  relations  compared  to  data  for  Series  D' 


Figure  7. 12.  Damage  prediction  relations  compared  to  data  for  Series  E' 
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Figure  7.13.  Damage  prediction  relations  compared  to  data  for  Series  F 
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Figure  7. 14.  Damage  prediction  relations  compared  to  data  for  Series  G’ 
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Figure  7.15.  Damage  prediction  relations  compared  to  data  for  Series  D'  including 
damage  initiation  adjustment 


Figure  7.16.  Damage  prediction  relations  compared  to  data  for  Series  E’  including 
damage  initiation  adjustment 
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Figure  7. 17.  Damage  prediction  relations  compared  to  data  for  Series  F’  including 
damage  initiation  adjustment 


Figure  7.18.  Damage  prediction  relations  compared  to  data  for  Series  G'  including 
damage  initiation  adjustment 
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Chapter  8 

SUMMARY  AND  CONCLUSIONS 

This  thesis  addresses  depth-limited  breaking  wave  damage  on  rubble  mound 
breakwaters.  Chapters  1-3  are  introductory.  Chapter  1  introduces  the  topic  of  rubble 
mound  structures  and  their  failure  modes.  It  is  noted  that  deterioration  models  are 
required  for  risk  and  reliability  analyses  as  part  of  rehabilitation  studies  for  these 
structures. 

Chapter  2  discusses  a  new  experiment  conducted  to  investigate  incipient 
motion  of  stone  armor  units  in  depth-limited  breaking  waves.  Observations  were  made 
of  incipient  motion  on  both  stone  and  sphere  armor  layers.  Velocity  measurements  were 
made  during  these  experiments  both  inside  and  just  outside  the  armor  layer  using  a  laser 
Doppler  velocimeter.  It  is  shown  that  the  maximum  vertical  convective  acceleration 
across  the  armor  layer  is  proportional  to  the  square  of  the  maximum  vertical  velocity. 
Equation  2.1 1,  an  incipient  motion  criterion  similar  to  the  Shields  criterion  for  sediment 
motion,  is  derived  based  on  Morison  forcing  for  the  dominant  mode  of  motion:  vertical 
lift  under  the  steep  breaking-wave  face.  Equation  2.1 1  is  restated  in  terms  of  the 
stability  number  as  Equation  2.12  as  follows 


144 


The  experimental  observations  and  measurements  are  used  to  validate  this  incipient 
motion  criterion  for  breakwater  armor. 

In  Chapter  3,  historical  damage  experimental  techniques,  measurements, 
and  models  are  described  for  rubble  mound  breakwaters.  It  is  noted  that  there  are 
inconsistencies  in  both  historical  and  recent  modeling  techniques  which  have  a  large 
effect  on  the  conclusions  drawn.  In  particular,  different  methods  of  computing  eroded 
area  or  volume  produce  results  that  differ  by  nearly  a  factor  of  2.  Several  authors  have 
made  note  of  the  fact  that  irregular  wave  experiments  may  not  produce  an  equilibrium 
profile  prior  to  failure.  This  conclusion  was  verified  within  this  study.  Many 
breakwater  stability  experiments  have  been  conducted;  but  nearly  all  tests  were 
conducted  with  nonbreaking  waves,  started  with  an  undamaged  structure,  and 
progressed  to  relatively  low  levels  of  damage. 


Chapters  4-7  discuss  a  new  experiment  to  investigate  damage  development 
on  a  breakwater  exposed  to  a  series  of  depth  limited  storms  of  varying  wave  energy  and 
in  varying  water  levels.  It  is  noted  that  this  study  was  motivated  by  the  fact  that  no 
general  irregular-breaking-wave  damage  progression  experiments  have  been  conducted. 
Therefore,  no  generally  applicable  relations  exist  that  would  permit  prediction  of 
damage  progression  in  breaking  waves.  A  methodology  is  developed  for  measuring 
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damage  in  a  small-scale  physical  model.  The  methodology  uses  damage  defined  by 
Equation  3.1  as  the  eroded  area  normalized  by  the  square  of  the  nominal  stone  diameter 


In  addition,  three  profile  descriptors  of  engineering  interest  are  suggested  as 
follows.  The  minimum  cover  depth  (Equation  5.1),  maximum  eroded  depth  (Equation 
5.2),  and  m^iT^irnum  eroded  length  (Equation  5.5)  along  the  profile  are  normalized  as 


T  -  =25 

~  D^'  e 

Data  were  obtained  from  an  automated  profiler  on  a  traditional  rubble  mound 
breakwater  section.  The  profiler  is  shown  to  provide  improved  damage  measurements, 
as  the  cross-shore  spatial  sampling  interval  was  less  than  1  mm,  which  is  considerably 
smaller  than  traditional  sounding  measurements.  In  addition,  the  profiler  was  automated 
making  profiling  very  efficient.  Because  profiling  did  not  require  draining  the  wave 
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flume  between  tests,  more  tests  could  be  accomplished  and  with  more  frequent 
profiling. 


In  the  experiment,  damage  development  was  measured  on  two  identical 
breakwaters  with  seaward  slopes  of  1:2  exposed  to  depth-limited  breaking  waves  on  a 
1:20  beach  slope.  The  armor  was  limited  to  a  traditional  two-layer  thickness  of 
uniformly  sized  stone  and  a  widely  graded  riprap.  The  waves  were  ran  as  sequences  of 
storms  with  varying  irregular  wave  conditions  and  water  levels.  Seven  unique  long¬ 
term  damage  test  series  were  conducted.  The  first.  Series  A',  was  a  long  duration  test 
with  varying  water  level  and  wave  height.  In  this  first  series,  the  water  level  was 
increased  midway  through  the  test.  The  wave  heights  were  increased  in  five  increments. 
Damage  was  allowed  to  reach  an  observed  apparent  equilibrium  level  before  the  wave 
height  or  water  level  was  changed.  The  second  and  third  series  (B'  and  C)  were 
designed  to  yield  insight  into  damage  progression  for  variations  in  storm  ordering. 
During  Series  B\  the  water  level  was  increased  midway  through  the  test,  while  for 
Series  C  the  water  level  was  decreased  midway  through.  For  Series  B'  and  C',  wave 
heights  were  increased  in  four  increments,  each  after  about  4,200  waves  or 
approximately  8  to  10  hr  of  prototype  storm.  Series  D’  and  E'  were  similar  to  Series  B' 
except  with  different  peak  wave  periods.  Series  F'  and  G’  were  similar  to  Series  B'  and 
D'  except  that  riprap  was  used  as  armoring. 

It  was  observed  during  the  series  that  the  damage  never  really  stabilized  or 
reached  an  equilibrium.  Also,  Series  5'  and  C'  indicated  that  the  final  damage  was 
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similar  for  both  series  consisting  of  different  sequences  of  storms  of  similar  cumulative 
wave  action.  The  damage  data  showed  partial  healing  during  damage  progression, 
where  eroded  pockets  from  displaced  stones  were  filled  in  by  minor  shifting  of 
surrounding  stones. 

The  probability  density  functions  of  damage,  eroded  depth,  and  cover  depth 
were  analyzed  in  Chapter  5.  For  each  series,  at  each  point  in  time,  the  probability 
density  functions  of  normalized  values,  given  by  Equation  5.6  as 


5*=^-^; 


are  shown  to  have  ranges  given  by  Equation  5.10  as 

-2.7<5*<3  ;  -2.7  <£'*<2.7  ;  -2.7  <  C*<  2.8 

All  normalized  variables  are  in  the  range  from  -3  to  3.  It  is  shown  that  series  with  lower 
overall  damage  levels  had  negatively  skewed  distributions  with  fewer  points  in  the 
negative  tail  in  the  range  -3  <  5*  <  -2  due  to  the  fact  that  there  were  few  measurements 
at  very  low  damage  in  the  range  5=0- 1. 


Equation  6.1  was  obtained  through  statistical  analysis  of  the  profile  data  to 


describe  the  alongshore  variability  of  damage  as 


This  equation,  along  with  Equation  5.10,  can  be  used  to  predict  the  range  of  damage  for 
a  given  level  of  mean  damage.  These  relations  indicated  that  alongshore  damage 
variability  was  greater  than  previously  suspected  and  could  be  significant  in  a  life-cycle 
cost  analysis. 

The  damaged  profile  shape,  as  characterized  by  the  mean  eroded  depth  E 
and  mean  cover  depth  C,  was  shown  to  be  well  described  in  terms  of  the  mean  damage 
by  Equations  6.2  and  6.5 

E  =  0.46 

C  =  Q  -  0.15 

where  Q  is  the  mean  cover  depth  when  S=0.  The  mean  eroded  length  was  shown  to 
be  described  in  terms  of  the  mean  damage  by  Equation  6.4  as 

I  -  4.4^5 

The  variability  of  eroded  depth  and  cover  depth  was  shown  to  be  described  by 
Equations  6.3  and  6.6,  respectively,  as 

o^  =  0.26  -  0.00007(5-7.8)^ 

+  0-098  -  0.002(5  -  If 
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where  Oco  is  the  value  of  Oc  when  S=  0.  It  is  noted  that  Equations  6.1,  6.3,  and  6.6 
underpredict  variability  of  erosion  somewhat,  but  more  data  are  required  to  better 
quantify  these  relations.  Equations  5.10  and  6.1  through  6.6  define  the  damage  profile 
shape  and  alongshore  variability  as  a  function  of  mean  damage  without  regard  to  wave 
conditions. 

The  temporal  variation  of  mean  damage  is  given  by  Equations  6.15  and  6.16 

S(t)  =  S(t„)  +  0.025 -  tf')  for 
S(t)  =  S(t„)  +  -  O  for 

for  time  domain  and  frequency  domain  wave  statistics,  respectively.  These  formulas 
can  be  used  for  predicting  mean  damage  S,  where  the  wave  height  and  period  vary  in 
steps.  The  formulas  fit  the  data  well  for  Series  A',  B\  and  C',  with  systematic  variations 
of  water  depth  and  wave  height.  The  equations  are  shown  to  fit  the  trends  of  damage 
progression  well  for  Series  D’,  E',  F,  and  G',  where  the  peak  wave  period  and  stone 
gradation  were  varied.  Damage  initiation,  where  only  1  or  2  stones  moved,  was 
consistently  underpredicted  by  more  than  a  standard  deviation  for  Series  D'  through  G'. 
This  appears  to  be  due  to  the  variability  in  damage  initiation.  It  is  shown  that  the 
prediction  is  significantly  improved  if  the  damage  progression  is  begun  immediately 
after  the  initial  profile  adjustment  or  1000  waves.  The  initial  profile  adjustment  may 
need  to  be  accounted  for  in  test  series  with  relatively  small  cumulative  damage.  Also, 
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the  empirical  parameters  in  these  equations  may  need  ,  to  be  varied  somewhat  with  wave 
period  and  armor  gradation. 

The  relations  given  in  this  thesis  for  eroded  area  and  damaged  profile 
prediction  can  be  used  to  predict  the  future  performance  of  a  new  or  damaged  structure. 
They  can  be  used  to  predict  damage  for  a  single  design  event  or  in  a  simulation  to 
predict  the  life-cycle  costs.  Further,  because  the  variability  in  damage  is  given,  these 
equations  can  be  used  to  predict  the  reliability  of  a  given  design  or  of  a  damaged 
structure.  The  equations  described  above  are  valid  for  the  range  of  wave,  water  level, 
and  structure  conditions  within  this  experiment.  Employment  of  a  critical  stability 
number  for  measurable  damage  may  be  necessary  to  reduce  the  number  of  storms 
required  in  a  life-cycle  analysis.  These  relations  should  be  conservative  because  the 
structures  were  relatively  impermeable,  relatively  low  crested,  and  exposed  to  severe 
breaking  waves.  But  they  may  not  be  conservative  for  more  highly  grouped  wave 
conditions  and  other  structure  slopes.  Therefore,  care  should  be  taken  in  extending 
these  relations  for  conditions  not  tested  herein.  Additional  experiments  are  being 
conducted  to  evaluate  the  generality  of  these  equations. 
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APPENDIX  A 


SUMMARY  OF  WAVE  MEASUREMENTS 
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Figure  A.l.  Wave  generator  command  signals  for  six  wave  cases  of  initial  experiment 
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Depth  =  11.9  cm 


Figure  A.2.  Spectra  of  wave  generator  command  signals  for  Waves  1,  2,  and  3  for  1 1.9 
cm  toe  depth 


Depth  =  15.8  cm 


Figure  A.3.  Spectra  of  wave  generator  command  signals  for  Waves  4,  5,  and  6  for  15.7 
cm  toe  depth 


154 


Eto 


Time  (seconds) 

Figure  A.4.  Wave  1  measured  time  series  from  wave  gages 
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Figure  A.5.  Wave  2  measured  time  series  from  wave  gages 
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Figure  A.6.  Wave  3  measured  time  series  from  wave  gages 
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Figure  A.7.  Wave  4  measured  time  series  from  wave  gages 
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Figure  A.9.  Wave  6  measured  time  series  from  wave  gages 
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Figure  A.  1 1 .  Wave  2  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  12.  Wave  3  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  13.  Wave  4  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  14.  Wave  5  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  15.  Wave  6  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  16.  Wave  7  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  17.  Wave  8  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  18.  Wave  9  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  19.  Wave  10  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.20.  Wave  1 1  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.21.  Wave  12  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.23.  Wave  14  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.24.  Wave  15  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.26.  Wave  17  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.27.  Wave  18  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.28.  Wave  19  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.  30.  Wave  21  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.32.  Wave  23  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.33.  Wave  24  incident  and  reflected  spectra  from  wave  gages 


S^(f)  (cmVHz)  (cmVHz) 


Figure  A.34.  Wave  25  incident  and  reflected  spectra  from  wave  gages 
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Figure  A.35.  Wave  26  incident  and  reflected  spectra  from  wave  gages 
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APPENDIX  B 


SUMMARY  OF  DAMAGE  AND  PROFILE  MEASUREMENTS 
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